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CYCLOTOMIC NUMBERS AND A CONJECTURE OF SNAPPER 


S. A. KATRE 


Department af Mathematics, S. P. College 
Pune 41} 030 


(Received 23 May 1988) 


The purpose of this note is to prove a conjecwre of Snapper [J. Algebra 97 
(1985), p. 277] by relating the numbers c;; defined by him with cyclotomic 
numbers and using elementary properties of Jacobi sums. 


1. INTRODUCTION 
Let F, be a finite field of g elements (9 not necessarily a prime) and let y bea 
generator of the cyclic group Fr. Lete > 1 be a divisor of g—1 and let © be a pri- 


mitive (complex) eth root of unity. Let % be the character on F, defined by X (y) = ¢ 
and % (0) = 0. For i, 7 modulo e, the e* cyclotomic numbers 4;, and the e? Jacobi 
sums J(i, /) of order e are defined by 


Ay cara. © F,\%(v) = oh (ytl) = 6} ie 
and 
Jui,jj)= 2% Xf (v) XI (v+1). rOR2) 
vera 
Let g = 1 + efand let Hy be the unique subgroup of ts of order f. In other 


words Hy is the subgroup generated by y°. Snapper® defines an eX e matrix C. = (ci) 
of nonnegative integers by 


oy, = Card. Hy O (x* +. 71.4 Ptkss) 
and conjectures that for fixed e, 
Cj > or asg >. ...(1.4) 


(See conjecture 8.1 in Snapper’). 


In section 2, we state some known properties of J (i,j) indicating their proofs. 
In section 3, we obtain the asymptotic behaviour of the cyclotomic numbers 4;,;. In 
section 4, we connect the numbers ci, with Aj; and prove Snapper’s conjecture. In 
section 5, we show that the ci, and the Aj, of order e are positive forg > e* (e > 1), 
a result of interest in connection with Section 8 in Snapper’. 
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2. ELEMENTARY PROPERTIES OF JACOBI SUMS 
A relation between Aj, and J (i, j) is given in 


Lemma 1— = & (-(4!+6)) J (i, j) = e? Aap. 
3 


The proof of this is similar to the proof of a particular case of this considered 
in, section 1 of Parnami et al.?. 


As in section 3 of Snapper’®; define u to be the unique integer such that —1€ y“Hy 


and 0 qu < e—1. Clearly X-) = @. By Proposition 3.2 of Snapper’, if Char. 


F, = 2,u = 0, If Char. F, > 2, then if fis evenu = 0; if f is odd, e is even and 
u = e/2. u may also be considered modulo e according to the context. 


Lemma 2— J (0,0) = q—2. For i, j $ 0 (mod e), J (i, 0) = — C™, 
J (0,j) = J(i, -i) = -1. 
Proor: J(0,0)= Zs 1 = q-2 
vEF,—{0, —1} 
J(7i,0)= & x0) = — x6 (-1) = — Gu 
v~-1 
JOj= 2 P+ I =-vV()==1. 
v #0 


J (i, —i) = = Xv) X#( + 1) = 


Paes, 
aorta te (a) 
= = x)= — x1) = -1. 
v1 


Lemma 3— J (i,j) J (i, j) = q for 
i, j,i + 7 0 (mod e). 


The proof of this lemma is similar to that of Lemma | of Parnami et al,*. 


3. THe Asymptotic BEHAVIOUR OF CYCLOTOMIC NuMBERS 


Let Xi, = VE F, |X) =0,X(v + 1) = &}, for i, J modulo e. Then 

At, = Card. X1;. Clearly the Xj, are pairwise disjoint sets and U X;, = F, — {0-I}. 
LAP 

Hence = & 4); = q—2, The cyclotomic numbers of order e are e* is number and it is 

rs 

expected that the elements of F, — {0, — 1} be almost equally distributed among the 

sets Xi, atleast for large g. This is confirmed by our following 
Theorem 1— The cyclotomic numbers of order e ar 

q~ co. (In particular the cyclotomic numbers tend to co along 

PRrooF : We have by Lemma 1, 


€ asymptotic to q/e? as 
with q.)} 
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e-1 
e* Aas = J (0,0) + & J (i, 0) C-o8 
{1 


+ 5° J(0, j) to 
j= 1 


e-1 
=F aca Theo) Lore ob. K AEM: 
where 
K= >> a (i, J) C7 (aitbs) 
i, j, i+j SE O (mod e) 
By Lemma 2, 
e? Ag = q—2 — > C-(a-u)t z aad) 
i=1 jal 
e=} 
=e itale-0)t 4 K 
i=1 
mg —2 —-(¢(a—u) —1) — (¢(d)—1) 
=(4(G—0) —1) + KR 
where 
peerenin | 4, 
< (@) fe if ea: } ---(3.2) 
Thus 
e? Aap = 9 + 1 — € (a—u) — € (b) — € (a—b) + K. ra oe) 
Here, 
|K|I< > lJ@ ij) | = (e-lI(-2)Vv¢a — .-- (3.4) 


i,j, i+j = = 0 (mod e) 


Dividing (3.3) by q and letting g — oo, we see that e? Aw/q > 1, i.€. Aas is asympto- 
tic to g/e? as required. 


4. PRoor OF SNAPPER’S CONJECTURE 


We first connect the c;, defined by Snapper (see (1.3)) with the cyclotomic num- 
bers A;,; in the following : 


Lemma 4— Ajj = Ctsusj3 Cty = Aisus;- 


Proor: A;, = Card. {re Pi lve y'Hyv+ti1€ y’ H; 


(equation continued on p. 102) 
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= Card. Wve Ee |v y-! © Ay, v yn! + yt yt) Ay }, 
= Card. ve F* lvy-€ Ay, vyfe — y-§ + WP Hy \, 


= Card. 1" E F* | veE Hy;, y«e— y! + ya! Hy \. 


= Card. Hy 0 (— y=! + 17 **! Hy), 
= Card. Hy 1 (7~? + y*-/ Hy) by Lemma 1.1 of Snapper’. 
= C_jt-j = City = Ci+u,; (by Theorem 3.1 of Snapper’). 
Hence 
Cty = Algyy; ASU = —u (mod e). 
From Lemma 4 and our Theorem | we have: 
Theorem 2— The entries c;; of the matrix C. are asymptotic to gle? as gq > ©¢. 
In particular, we have, for fixed e, 
Ct; > ©O aS g—>oo ...(4.1) 


proving the conjecture of Snapper. 


5. A REMARK 


In Proposition 8.1 of Snapper’ it is shown that certain Fermat-type equations 
have no nontrivial solutions in F, (q a prime) or in integers, provided ci; = 0. As we 
have shown that for fixed e, Ci; * °° aS gq > oo, it may be interesting to see after 
what stage the c;, (and so also-the cyclotomic numbers) become positive. 


From section 5 of Snapper*, we see that the ci; and the Aj, are positive for e = 1 
when g > 2, and for e = 2 when q > 5. From (3.3) and (3.4) we get 


e? Aab > 9 + 1 — € (a—u) — (b) — « (a—b) — (e—1) (e—2) V @. 
certo) 
Here gq is not necessarily a prime. 


(Compare this with the result of Dickson’ [45], viz. for q = P, a prime, and e an 
odd prime such that p = 1 (mod e), 


e* Aoo > p — 3e +1 — (e~1) (e—2) Vp) ...(5.2) 


From (5.1) we get, for all i,j (mod e), 


e* As 2 g'— (e~1) (e—2) Vg — (Be—I). (5.3) 
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Now, for x, b, c > 0, we have, 
x* — bx — c > 0 provided x? > b? + 2c, ---(5.4) 
Hence the right hand side of (5.3) is positive provided 
q > (e—1)? (e—2)2 + 2(3e—1). BA RP) 
Therefore we have, 


Proposition 1\— The cyclotomic numbers Aj, and the numbers c;; of order e are 
positive for 


q > (e—1)* (e—2)? + 2(3e—1). 


The right-hand side of (5.5) may be written as e — e (3e—2) (2e—3) + 2. Hence for 
e > 2, the right-hand side of (5.5) is less than e*. Thus we get a slightly weaker but 
interesting result viz. 


Proposition 2— The Aj; and the c;, of order e are positive for g > e* (e > 1). 


A result sharper than Propositions 1 and 2 but more awkward looking, which 
follows from (5.1), is 


b+ Vbt+ 54) 


Proposition 3—The Aj; and c;, of order e are positive when g > ( 5 


where b = (e—1) (e—2) and c = 3e —1 or 2e—1 according as u = 0 or e/2. 


The bounds in these propositions can be sharpened to some extent for particular 
values of e using cyclotomy. 
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STABILITY IN MAMMILARY COMPARTMENTAL SYSTEMS 


H, EL-Owaipy, A.A. AMMAR AND O.A. ELLEITHY 


Department of Mathematics, Faculty of Science 
Beni-Sweif, Egypt 


(Received 26 July 1988) 


This paper is concerned with the stability of mammilary compartmental 
system with constant transfer rates and no ingestion of material to any com- 
partment from the outside environment. 


1. INTRODUCTION 


Compartmental analysis has recently been one of the most widely used mathe- 
matical techniques in biological modeling (see for example Jaquez‘). In this technique 
a system is divided into “compartments” with laws defining rates of exchange between 
them, usually exchange of some material. In physiological models*’® the compartments 
could be suborgans of an organ in the body, whereas in ecological model’, the com- 
partments could be populations. 


In the usual representation for a compartmental system, a box (ora point) 
denotes a compartment, and an arrow indicates the transfer of material into or out of 
a compartment (see Fig. 1), There also may be inputs from the outside environment 
into one or more compartments (vertical arrows pointing into the tops of boxes) and 


¥ 


Fic, 1, 
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there can be excretion of material from some of the compartments to the outside 
environment (vertical arrows pointing out of the bottoms of boxes). Should there be 
no exchange of material to the outside environment, the compartmental systems are 
referred to as closed, otherwise, it is said to be an open system. Realistically, many 
compartmental systems are open, for some material are lost due excretion, metabolism, 
etc. 


The goal of this article is to investigate the approach to equilibrium point 
(stability) in mammilary compartmental systems (this will be done in Section 3). A 
mammilary compartmental system has a central compartment called “mother” which 
exchanges material with all the other, ‘daughter’? compartments, but there is direct 
exchange between any two of the daughters. Each cycle of length exactly two. Levine? 
considers an example (a nonlinear closed mammilary compartmental system) of a 
mammilary compartmental system. He considers the resources R as “mother” and 
n species as daughters with biomasses N,, N,, ..., NV, which lose biomass to the en- 
vironment at exponential rates a,, ..., a, respectively. Then, N,, ..., N, obey equations 
of the form 


N= — a N+ fi(K — Nine Dimsleng n 
=) 


Sins Sita > 0 


where /; are the ingestion functions from the source R into the compartment i, fj; and 
fiz denote partial derivatives of /; with respect to the first and second arguments, and 


K= 2 N, + R (i.e. the total biomass is constant). Daniel S. Levine shows that the 


local stability of a positive solution (if it exists) under the condition that the ingestion 
functions f; (i = 1, ..., 1) are grown slower than linearly with N;. He derives some pro- 
positions concerned with competitive systems. Finally he put the following question : 
“whether approach to equilibrium occurs for all mammilary systems or more parti- 
cularly for mammilary systems with structure described in his article ?”. 


2. THE MopDEL 


Let x1, ..., X, denote the biomass of compartments 1, ..., respectively (see for 
example Fig. 2) with central vertex (mother) compartment | which exchanges material 
with all the other compartments (daughters) such that there is no direct exchange of 
material between any of two daughters and we assume that there is no input of 
material from the outside environment. The biomasses x;, ..., X, (here, the variables) 
obey the equations : 

(n-1) 
&, = 01, %1 + fer by Xt4; 


(*) 


&, = a, X, — 5, X, 
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hy = A,X, — 5, Xs 
x, = Qn-1 xX; es b,-1 Xn 


} stants and eis sory bk are the 
wher ents a, b; are nonnegative con ; 
ere the transfer componen ’ : 
tes of change of biomass of the compa! tments |e cosy M respectively. This model has 
ta 


matrix of the form : 





| 41, b, 5b, bn -1 | 
. a —b, O 0 | 
ee ee BoP oe] 2.1) 
A= : 
Ae aoe: 
{ An-1 n—L j 
(n-1) (n-1) F 
where a1, = — 2% aj, if the system (*) is closed and a1, = — [a,, + Re ai] if the 


system (*) is open, a), is an excretion of material (biomass) from compartment a | 
(mother) to the outside environment. In the case of open mammilary compartmental 


systems, the trivial solution is locally asymptotically stable as shown in the following 
theorem : 


3. THEOREM 


Theorem— The trivial solution of mammilary compartmental system (*) is 


locally asymptotically stable if there exists an excretion of material from some of the 
compartments to the outside environment (that is if the system is open). 


Proor ; The eigenvalues of the matrix A (2.1) are real b 
is of length which does not exceed two 
(see H. Anderson’, p. 61). Such matrix 


ecause each of its cycle 
(in this case each cycle of Jength exactly two) 
must have nonnegative off diagonal elements, 
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then if the determinants of its principal minors alternate in sign, the eigenvalues of A 
have negative real parts’, but they are Originally real', then we can deduce that the 
eigenvalues of A are real and negative. Consequently, the trivial solution of the system 
(*) is locally asymptotically stable. To verify this we restrict ourselves (for simplicity) 
to the mammilary system with five compartments (see Fig. 2), its matrix has the form: 


hae Dy Oy be bAA) 


ia—b 0 0 0 | 
A= ja O0-b 0 mI pend Oc) 
| fae ot bs | 
Ea 8 aCe he 
where a,, = — Ba if the system is closed, led dj, d, ds, d,, d; denote the deter- 


minants of the principal minors of A respectively. If the system is closed the verifica- 
tion of the stability is difficult since d; = | A | = 0 (in this case: d= 0) de>, 
d; < 0,d, >0,d; =| A|=0). For this reason, we assume that there exists an 
excretion of material from some of the compartments to the outside environment, for 
obtaining d; < 0 to conclude that the trivial solution is locally asymptotically stable. 
Thus, the following procedure restrict to check the sign of the principal minors of A 
in (3.1) (under the assumption that a, 5; and ay, (i = 1,2, ...,”) are positive cons- 
tants, where do: is an excretion of material from compartment, i, to the outside 
environment) : 


(i) The case of an excretion a9, > O of material from the compartment | to the 
4 
outside environment (see Fig. 3), in this case a,; = —dy1 — Pah a;, One can show that 


the determinant of the principal minors of the matrix A are: 


3 ye 


2 a,, 5 


Fic. 3. A five vertices mammilary systems with all excretion in compartment 1. 


Ij 


4 
d, yd aaa [ao, + pede <0 


d, = [a, + a3 + a4 + ap,] db; > O 
d; = — [a3 + a, + a,] bd, b2 < 0 
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d, = [ay + ao:] 0; 5, 6; > O 
d,=|A|= — di 5, b, b, by < 0. 


(ii) In the case of an excretion @y2 of material from compartment 2 to the out- 
side environment (see Fig. 4). In this case, the matrix A has the form : 


[ au by be bs b, } 
—b,—a,, 0 0 0 | 
| & 1 02 
Ace | Ban bed Oh SOG Bs eo 
| a3 0 0 —b, UO | 
t a, 0 0 06—57) 


; + 


2 5 
Ao2 


Fic. 4. A five vertices mammilary system with excretion of material from compartment 2 
to the outside environment. 


4 
where a), = — 2 a;. One can show that the determinants of the principal minors of 
A are 
4 
{=] 
d, 


= [a, + a, + a] (b, + doz) + @1 ao > 0 


dy = — [as by + ay bi + (a, + as + a) age) b, < 0 
d, = [(a, + as) G2 + 44 b)] b, bg > O 
ds=|A|= 


— A, Ags by bs by < 0. 

(iii) In the case of an excretion as > 0 of material from the compartment 3 to 
the outside environment the determinant d:, d, of the first and the second principal 
minors are the same as in the case of closed system i.e. d, = a,, < 0, d,=(a,+a3+ ds) 
b, > Othen it remains to check the sign 


: of the determinant of the other principal 
minors. The matrix A in the case has the form : 
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f i Dee d, bs by / 
| a, —b, 0 0 0 | 
A= |a O0-b,-aq30 O | 
| a, 9-0 —b, 0 | 
iP 0 0 0 —bd, ; 
d; = — [a b, + da bz + a, ays + a3 Ogg + ay yg] 6; < 0 


d, = [a4 b, -}- as Qo3 + a, Qos] by bs > 0 
d; = | A | = a, Qo3 b, os ba < 0 
(iv) If there exists an excretion aos of material from compartment 4 to the out- 


side environment it suffices to examine the sign of d,, ds. The matrix A in this case 
has the form : 


ox bb, bs by | 
A= | as 0 —b, clei! 8, 
| 
| 
L 


j}a3 O O —b3—ao, 0 

OPscal 0! 70TEBel 

One can show: 
ds 

d; 


[a4 bs + QA, Ao4 + das aoa] b, b, s.0 


I 


\\ 


| A | = — a3 doy 0, bz bg < Oz 


(v) Finally, if there exists an excretion a), > 0 of material from compartment 5 
to the outside environment, the determinants di, d,, dz, dy are the same as in the closed 
case then it suffices to check the sign of the determinant of the last principal minor : 


d. ase 44 Qos D, b, b, <0. 


o 


Thus, we observe that in each of the previous case the sign of the determinants of the 
principal minors of the coefficient matrix A alternate in sign then the eigenvalues of 
A have negative real parts’, but the eigenvalues of A are originally real’, then it 
follows that the eigenvalues of A are real and negative proving the local stability of 
the trivial solution for mammilary system with five compartments, it is plausible 
to generate this result to the mammilary systems with n compartments. 


4. DISCUSSION 


In this paper we discussed the stability of mammilary compartmental system 
(with constant transfer rates) : 


x = Ax 
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with no ingestion of material to any compartment from the outside environment. It 
remains to study the stability of mammilary compartmental system with nonzero 
ingestion of material, b;, into any compartment, /, from the outside environment, the 
dynamics of this system with constant flow rate are based on a differential equation of 


the form : 
k= Axc+b 


where x and b arenx1 column vectors, A is nxn compartmental matrix. We leave 
this study to a future work. 
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The problem of determination of an optimal programme for augmentation of 
capacities of depots and the number of buses to be parked overnight at 
respective depots and their shipment from depots to the starting points of 
routes is considered. The objective is to minimize the total capital expendi- 
ture to be incurred in augmenting capacities of the depots plus the present 
value of the total cost of travel performed by buses between the depots and 
the starting points of the routes over a planning horizon. A method is 
developed to obtain the solution of this problem. Further, a software pack- 
age for the method is developed and tested on ICL 2900 Computer. 


1. INTRODUCIFION 


Problems relating to urban bus transportation system have been of interest for 
quite some time. 


In the past, these problems have been tackled by applying rudimentary techniques 
based on commonsense and experience. This has resulted at times into avoidable 
expenses. On account of this and growing complexity of problems relating to urban 
bus transportation system and increasingly larger sums of money involved, the need 
to employ analytical tools to deal with these problems has been lately felt. In recent 
times, analytical tools have been attempted to deal with them by several resear- 
chers*®. 

Buses are parked overnight at depots. Starting points of routes are generally 
different from depots and a bus has to travel from its depot to the starting point of 


its route before it can be engaged on regular service. Similarly, at the end of the 
service, the bus travels back to its depot. Distance traversed by a bus in going from 
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its depot to the starting point of its route and back pays no return and is known as 
dead mileage or dead-travelling. It is desirable to reduce dead-travelling. On the 
‘other hand, depots built in the past cannot provide parking facilities to the increased 
number of buses necessitated by an increase in the number Of city-bus travellers. 
Capacities of the depots have to be augmented to satisfy the enhanced demand of 
parking facilities. Capital expenditure to be incurred in creating capacity for a bus 
varies from depot to depot. There is an upper bound on the number of buses by which 
capacity at each depot can be augmented. The problem is to determine the number 
of buses by which capacities of different depots should by augmented and the number 
of buses to be parked overnight at them and the number of buses to be sent from each 
depot to the starting point of each route with an objective to minimize the total 
capital expenditure to be incurred in augmenting capacities of the depots plus the 
present value of the total cost associated with the dead-travelling of buses over the 
planning horizon. This problem is formulated as a capacitated transportation-type 
problem and is solved by the procedure given by Dantzig? after making some altera- 
tions. The solution procedure is illustrated through a numerical example and tested 
on ICL 2900 Computer. Utility and versatility and future extension of the present 
work are also discussed. 


2. FORMULATION OF THE PROBLEM 


Suppose that there are m existing depots, p sites for potential depots and n 
routes. Let the existing capacity measured in terms of buses at depot i (i = ], ..., 
-m-+ p) bea, Fora potential depot i (i = m+ 1, ..., m + p), the-existing capacity 
a; is equal to zero. Let 6, (j = 1, ..., ) be the current enhanced demand of bused on 
route j, c the cost in rupees of running a bus per kilometre, c; (n41) i = 1, ..., m + p) 
the expenditure in multiple of thousand rupees to be incurred in augmenting capacity 
for one bus at depoti, dj, (i= 1,...m+p;j7=1,..., n) the distance in kilo- 
metres between depoti and the starting point of route Js fitnaty G@ = 1, 3... m + p) 
the upper bound on augmentation of capacity of depot i, Nthe planning horizon in 
years, r the annual rate of interest, x; (i = 1, .... m + p) the number of buses to be 
parked overnight at depot i, xi, (/ = 1, ....m + p;j = 1,..., n) the number of buses 


to be sent from depot i to the Starting point of route j, x; (n+ 1)G=1,. 


<5 NM + 
the augmented capacity at depot i. * 


The objective is to minimize the total capital 
expenditure to be incurred in augmenting capacities of the depots plus the present 
value of the total cost associated with the dead-travelling of buses between the depots 
and the starting points of the routes over the planning horizon : 
formulation of the problem is as follows. . 
m+ p,j = 1,..., + 1) which minimize 


The mathematical 
Find integer x; and xi, > 0G — ],..., 


m+P n+l 


pm SS Ss bce (1) 


i=l j=l 
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where 
N 
Cha | 365 (2) dj, >. (1 + r[100)~é ] 1078 
t=1 
(fae do ce ge) 2 b> a) ee? 
subject to the constraints 
io Xt4j — Xt) =a (i = 1,...,.m + p) <52(3) 
jal 
m+P 
x, = by (j = 1,...,n) 
— ...(4) 
N 
x= Sy X14; (i = |, ...,m 4 p) real pe 
Yea) | 
and the upper bound restrictions 
Mi Gnai) S fina) Geslog; m+ Pp). ...(6) 


It may be noted that o,(i = 1,...,m +p;J=1,...,n) in multiple of 
thousand rupees is the present value of the cost associated with the dead-travelling of 
a bus between depot i and the starting point of route / over the planning horizon. 
Equation (5) are the consequence of the fact that the number of buses to be parked at 
depot i should be equal to the number of buses being sent from depot / to the starting 
points of the various routes. 


3. SOLUTION PROCEDURE 


The problem formulated above can obviously be solved using the branch and 
bound method of integer programming. But as this method involves much computa- 
tional work, an efficient method based on the procedure for solving the capacitated 
transportation problem is developed. 


Note that eqns. (5) are not active constraints of the formulated problem. So an 
optima! solution of the formulated problem can be obtained ignoring eqns. (5). The 
optimal solution of the formulated problem ignoring eqns. (5) would determine 
required values of all x;,’s. After the values of the x;,’sare determined, eqns. (5) 
are used to compute required values of all x;’s. As the formulated problem ignoring 
eqns. (5) possesses total unimodularity for its coefficient matrix and a;’s and b, ’s are 
integers, integrality of all x:,’s of its basic feasible solutions is guaranteed. Unimo- 
dularity and its implications are discussed by Bazaraa and Jarvis'. With the above 
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observation, the formulated problem ignoring eqns. (5) is a capacitated transportation- 
type problem which differs from the usual capacitated transportation problem con- 
sidered by Dantzig® in that the coefficients of x; (nm + 1) (i = J, ..., m™ + p) in con- 
straints (3) are —1 init, whereas the coefficients of all the variables appearing in the 
constraints are | in the usual capacitated transportation problem. To obtain the 
optimal solution of the formulated problem ignoring eqns. (5), work proceeds as 
follows. 

To overcome the difficulty in obtaining an initial basic feasible solution of the 


formulated problem ignoring eqns. (5), we introduce artificial integer variables 
X¢mypei); 2O (J = 1, ..., n) into eqns. (4) which then assume the form 


m+P+1 
x, = 6,(j = 1,..., ”). 7) 
i=l 
The upper bound restrictions on aritficial variables are 
X(mtp+t)j So (J = 1, ... nn). (8) 


Entry of artificial variables into optimal solution is prevented by associating a 
cost M with each of them where M is an arbitrarily large positive number. After this, 
the formulated problem ignoring eqns. (5) reduces to the following augmented 
problem which requires determining integer x1; > (i = 1,...,.m + p+ 1;j7=1,.., 
n+ 1); but noti = m +p + 1 and J = n + | simultaneously) that minimize 


m+p n+) 


a = a Ss Ci; Xij a M Si X(m4p+1)j (9) 


i=l jai ja) 


subject to the constraints (3) and (7) and upper bound restrictions (6) and (8). The 
augmented problem also possesses total unimodularity for its coefficient matrix and so 
integrality of all x, ’s of its basic feasible solutions is guaranteed. Solution of the aug- 
mented problem provides the soultion of the formulated problem ignoring eqns. (5). The 
tableau representation of the augmented problem is shown in Table I. In this Table, 
the left top corner of each cell (m + p + 1, j) in the row with the heading ‘Artificial 
variable” depicts the cost associated with the artificial variable x (m + Pipe 
while the left top corner of each cell (i, + 1)in the colu 
“Augmented capacity” 
for one bus at depot i. 


mn with the heading 
depicts the expenditure to be incurred in augmenting capacity 
: And the left top corner of each other cell (i, 7) in the tableau 
depicts the present value of the cost associated with the dead-travelling of a bus 
between depot i and the Starting point of route j/ over the planning horizon. The left 
middle space of each cell (i,m + 1) in the column with 
capacity” depicts the upper bound on augm 
the left middle space of all other cells depic 


; the heading ‘Augmented 
entation of the capacity of depot i. And 
ts anentry of co indicating that there is 
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TABLE | 
Augmented problem 

Starting points of Augmented Existing 
route 1 route 2 ae route n capacity capacity 

Xila +1 ) a; 

a ae ee ee se a ee ee ee ee ee 
Cj Ciz Cin Ci(m+1) 
Depot 1 i aC wes = fi ay 


1 1 <i) 











C21 Ca Cen €2 (n+1) 
Depot 2 o ~ 3 a fa (n41) a: 
1 1 1 —1 
C (m+n)? C (m4p)* C(mip)n C (m+n) (n+1) 
Depot (m+ p) 9 a vs f (m+n) (n+1) a mip 
1 1 1 —1 
Artificial M M M 
variable = - a = 
X (m+n+41)5 1 1 1 
Buses 
required bj; db b, Dd» 





no upper bound restriction on values of the variables corresponding to them. To 
obtain a; in the row with the heading ‘Depot i’, sum the products obtained by 
multiplying x;, with the entry in the left bottom corner of the corresponding cell 
across the row and to obtain 5, in the column with the subheading “route j”, sum the 
x1,’sacoss the column. And Z’ is obtained by summing the products obtained by 
multiplying x;; with the entry in the left top corner of the corresponding cell all over 
the tablaeu. 


An initial basic feasible solution for the augmented problem is found almost in 
the same way as for the standard cost minimizing transportation problem. Several 
methods for obtaining an initial basic feasible solution for the standard cost minimiz- 
ing transportation problein are discussed by Hadley*. Any of these methods after 
some modification can be used to obtain an initial basic feasible solution for the 
augmented problem. In the case of the standard cost minimizing transportation 
problem, all the methods for determining an initial basic feasible solution assign a 
nonnegative value to a variable and at the same time satisfy either a row or a column 
constraint at each step. But in the case of the augmented problem, a nonnegative 
value is assigned to a variable without forcing it to exceed its upper bound and at the 
same time satisfying either a row or acolumn constraint. If the value of the variable 
is limited by a row or a column constraint, the variable is considered basic. If, on the 
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other hand, the value of the variable is limited by its upper bound, it is considered 
nonbasic. However, if the value of the variable is limited by its upper bound and 
also by a row or a column constraint, it is considered basic. It is worth noting that 
x;; multiplied by the entry in the left bottom corner of the corresponding cell (i,j) is 
the amount used of resource a; and x;, alone is the amount satisfied of requirement 
b, in the case of the augmented problem, whereas x;,; is the amount used of resource 
a, and satisfied of requirement b, in the case of the standard cost minimizing trans- 
portation problem. The value of a basic variable is entered in the right bottom 
corner of the associated cell after enclosing it in a circle. Alsothe value of a non- 
basic variable different from zero is entered in the right bottom corner of the associa- 
ted cell but after making a bar over it. 


After a basic feasible solution has been obtained, it is tested for optimality. To 
do this, the relative cost coefficients ¢;,’s corresponding to the nonbasic cells are 
computed and their values are entered in the right middle spaces of these cells. To 
compute the ¢, ’s we proceed as follows. The(m-+ p)u;’s and nj’s are computed 
from the following three different forms of equations 


(corresponding to the basic cells in } 
| rows with the heading of depots | 
‘and columns with the subheading | 
| of routes 
| corresponding to the basic cells | 
lin the column with the heading | 
| of augmented capacity | 
yj ot | corresponding to the basic cells | 
|in the row with the heading of | 
| artificial variable | 


uM t+ VY, = C1; 


uy = —c,(n + 1) ...(10) 


Once the (m + p) ui’s and n v,’s are known, @j;’s are computed from the following 
three different forms of equations 


{ corresponding to the nonbasic } 
Cy =O, — u-V, | cells in rows with the heading of | 
|depots and columns with the! 
| subheading of routes 


| corresponding to the nonbasic 
| cells in thecolumn with tbe head- 
| ing of augmented capacity 


CH (n4-1) = Cina) tut 


ENTS 


ae ae ee ee 


| . 
| Corresponing to the nonbasic cells | 


lin the row with the heading of | 
| artificial variable | 


Ci (m4p4.1)) = M— vy, 
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If all the relative cost coefficients corresponding to the nonbasic variables at zero 
level are nonnegative and at their upper bound are nonpositive, the basic feasible 
solution is optimal. On the other hand, if at least one of the relative cost coefficients 
corresponding to the nonbasic variables at zero level is negative or at their upper 
bounds is positive, the basic feasible solution is not optimal. Then the nonbasic variable 
with which is associated the greatest numerical value of the relative cost coefficient 
among all the nonbasic variables at zero level with negative values or at their upper 
bounds with positive values of the relative cost coefficients is selected to enter the 
basis of the new feasible solution. Further, ifthe nonbasic variable selected to enter 
the basis is at zero level, 6 > O is added to its value and values of the basic variables 
are so adjusted that all the row and column constraints are satisfied. Then a value 
as large as possible is assigned to @ so that the value of no variable exceeds its upper 
bound. This yields a new basic feasible solution. And if the nonbasic variable 
selected to enter the basis is at its upper bound, @ > 0 is subtracted from its value and 
then the same procedure is followed as in the case of the nonbasic variable at zero 
level selected to enter the basis. Irrespective of what type of nonbasic variable enters 
the basis, the leaving basic variable is chosen as the one which reaches eithera zero 
value or a value equal to its upper bound. To ensure unique selection of leaving 
variable, that variable among the potential leaving variables is chosen as the leaving 
variable for which row index i is minimun. In case, there are two or more potential 
leaving variables with the smallest row index, then the one with the smallest column 
index jis chosen as the leaving variable. The new basic feasible solution thus 
obtained is again tested for optimality and the procedure is repeated till the optimal 
solution is obtained. Finally, it may be mentioned that if the optimal solution of the 
augmented problem has one or more artificial variables with positive values, then the 
formulated problem has no feasible solution. 


4. A NUMERICAL EXAMPLE 


Now we shall apply the above procedure to obtain the optimal solution of a 
numerical problem which is obtained by taking m = 2,p=1,n=4,c = 4, 
d,, = 5, di, = 7, dy = 3, dg = 4, dy = 8, doz = 9, dhs = 10, dg = 4, ds, = 6, 
dz. = 2, dss = 8, d34 = 6, N = 10,r = 10, and assigning numerical values to all 
other quantities in the problem formulated above in section 2. The tableau repre- 
sentation of the augmented problem associated with the numerical problem is shown 
in Table II. For this augmented problem, the objective function which we seek to 
minimize is 


(90x1, + 126 x5 + 54 x,3 + 72 X14 + 630 X45 } 
Z' = | + 144 x,, + 161 x2. + 179 x23 + 72 X24 + 550 X25 | 

~ 4+ 108 xs, + 36 X22 + 144 x5, + 108 xge + 610 x25 7 
Lae M (x4, + X42 + X43 + X44) J 


(12) 
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TABLE II 


Augmented problem Associated with numerical problem 


i a ne 





Starting points of Augmented Existing 
route | route 2 route 3 route 4 capacity capacity 
xi (441) a; 
90 126 54 72 630 
Depot 1 nts 7 = oe 35 35 
1 30 15 1 1 —1 
144 161 179 72 550 
Depot 2 eo 00 eo Se 50 30 
1 1 20 1 10 1 —1 
ee ges et ny to Se ed Pee ee oe Reread or aera eee 
108 36 144 108 610 
Depot (2 + 1) eo oo 90 ad 20 0 
1 eo 1 1 =] 
Artificial M M M M 
variable oo 20 °° 0 
X (24141)5 1 1 1 30 1 635 
Ad he” eee We ER tage OP 
Buses 30 25 40 35 
required b; 
SSNS 
TABLE IIf 


Final tableau providing optimal basic feasible solution 


Starting points of 





Augmented Existing 
route 1 route 2 route 3 route 4 capacity capacity 
X; (441) ai 
“ey |, 290 > TRG Ip Apsaah © rg ae ee 
Depot 1 end 2° Pag se 2° 680 35 35 
1 1 1 40 1 —{ <5 
=i) ==> G4 > genial “ahinaaegaiyg ai Rosey anne 
Depot 2 2 aad oT AS ~ 50 30 
1 30 i 1 inva —-1 40 
086 ae 
Depot (2 + 1) — 69 Se vate kK) TE aa 20 -65 0 
1 1 20 1 1 -1 20 
Artificial M M M M 
Variable °° M-694 © M-711 © 44.684 © M-622 
X (o4n41)5 1 1 1 1 / 


ee 
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Following the column-minima method after some modification as explained in 
section 3, an initial basic feasible solution of the augmented problem associated with 
the numerical problem is obtained and is shown in Table II. Skipping the inter- 
mediate steps, the final tableau providing the optimal basic feasible solution is shown 
in Table III. In Table III, values of all the relative cost coefficients corresponding to 
the nonbasic variables at zero level are nonnegative and at their upper bounds are 
nonposilive indicating that the basic feasible solution is optimal. 


According to the optimal solution, capacities of depots 1 and 2 should be 
augmented to have 5 and 40 more buses respectively. A new depot should be con- 
structed at site (2 4- 1) to provide parking facility for 20 buses. All the 40 buses 
should be sent from depot | to the starting point of route 3. 30, 5 and 35 buses 
should be sent from the newly constructed depot (2 + 1) to the starting point of 
route 2. The minimum capital expenditure to be incurred in augmenting capacities 
of the depots plus the present value of the total cost associated with the dead-travell- 
ing of buses between the depots and the starting points of the routes over the planning 
horizon is 47875 thousand rupees. 


5. RESULTS AND DiSCUSSIONS 


Computer programmes for the solution procedure were developed and tested on 
ICL 2900 Computer. The solution procedure based on capacitated transportation 
method was found to be computationally more efficient than that based on branch and 
bound method of integer programming with respect to Central Processing Unit execu- 
tion time and number of iterations. 


6. VERSATILITY AND APPLICATIONS OF THE MODEL 


The mode) considered is quite versatile and has applications in areas other than 
urban bus transportation system. This can be seen if the terms ‘depot’ and ‘starting 
point of route’ used in the present work are stretched to include a supply point and 
a demand point respectively. With this generalization, the model considered can be 
employed to find optimal solutions to similar problems arising in areas of food grain 
movement from godowns to distribution centres and of crude petroleum supply from 
production sites to refineries. 


7. Scope or Futurs Work 


One of the areas for future work is the development of a technique superior to 
the existing branch and bound technique of integer programming and the one develo- 
ped above for solving the problem considered in this work. This may be accompli- 
shed by exploiting the fact that the number of depots is quite small compared to the 
number routes in the problem considered in this work. The other area for future 
work is to solve the problem considered in in this work by introducing expressions 
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for costs of augmentation of capacities of depots which could take care of an initial 
fixed expenditure associated with the augmentation of capacity of a depot. 
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In this paper we study the existence, uniqueness, error estimations and con- 
tinuous dependence of solutions of a class of nonlinear higher order differential 
equations with the given initial conditions. Our approach is based on con- 
verting the equations into equivalent integral equations and the application 
of the Wazewski’s general method of successive approximations. 


l1. INTRODUCTION 


Let n > 2 be an integer and let p; (t), 0 < i <n, be positive continuous fun- 
ctions on J = [fo, a], fo > 0 and a > O is finite but can be arbitrarily large. We define 
the differential operators Z;,0 < i qn, by 


x (t) 1 4d 
Lbx(t)= — ~Ux(t) = ——sL-,x ji 
0 X (t) i (t) mie) Tei x, 1 Sign. 
In this paper we consider the nonlinear differential equation of the form 
CL, x(t) = F(t, Lo x(t), L; x(t), ..., La-1 x (t)) -»+(1) 


with the initial conditions 
fa-4, (fo) = Ci. ed Pe ta) 


where F: J X R" > Ris acontinuous function, C;-; are given constants and R de- 
notes the set of real numbers. 


There are many papers written on the various special forms of (1) — (2) from 
different points of view, for example, see’*-*!'!-!® and some of the references given 
there. In particular, Fink and Kusano* and Trench’’ have considered eqn. (1) and 
obtained conditions which imply that equation (1) has a solution x which behaves for 
large positive f like a given solution z of the unperturbed equation L, z = 0. The main 
results in®?! are established by using the well known Schauder-Tychonoff fixed point 
theorem. The object of this paper is to establish some results on existence, uniqueness, 
error estimations of solutions of (1) — (2) and also continuous dependence of solutions 
on the right side of equation (1). Here we study the problem (1) — (2) by converting 
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it into an equivalent integral equation and by using the general method of successive 
approximations based on the idea used by Wazewski'* (see, also Kwapisz and Turo’, 

~ Pachpatte?’). Our formulation of the more general problem (1) — (2) is motivated in 
part by the studies of various types of equations by the authors’®?*181516 and the 
recent papers of Fink and Kusano® and Trench". 


2. STATEMENT OF RESULTS 
We first convert the problem (1) — (2) into an equivalent integral equation. We 
say that x is a solution of (1) — (2) if Zo x, ..., L, x exist and satisfy (1) — (2) on J. 
We shall be interested in solutions x which are continuous in J together with Lp x, ..., 
L, x. The set of all such solutions will be denoted by C* (/). If a function x € Ges) 
is a solution of (1) — (2), then for the function » continuous on J and defined by the 
formula y (t) = L, x (t), we have 


La-1 X (t) = qn (t) + 1,» (t) 
Li-s x(t) = g,-1 (t) + Jy-1 Y(t) 
Li x(t) = q(t) + ny (t) 


Lox(t)=qn(t)+hy(t) .--(3) 
where 
Yn (t) a Cos 
‘ 
Qn-1 (t) — C,-: af Cee J Pn-1 (t,-1) dt, 
. ‘0 
ga(t)=C,+C, J Pz (ts) dtz + ... 
0 
t 's oak 
+ C,-) J P2 (ta) J vee J Pa-\ (tn) dt,-1 wes dt, 
"9 fo ‘o 
M(t)= Cot Cy f pi (t,) dt, + ... 
o 
‘ a Est 
T Cn—1 J Pr (t)) J “igsh ie Paes (Ean) Aber iiccentlty 
and : : ‘0 


I, y (t) = i Pa (S) y (s) ds 


0 


t 


L-ay(t)= J Pn-1 (tn-1) ue Pn (5) y (Ss) ds dty_, 
0 


0 


(equation continued on p. 122) 
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Tey) = F welt) Fone T* pacr ead Via (8) 9 (8) ds dtpyn dt 
0 ; 


0 0 0 


t 
1 


‘ q P= 
y(t) = J P (ti) J er A Bag ea) Pp, (s) » (s) ds dt,_1...dt,. 
0 


7 
0 0 "9 
Consequently the function y fulfils the equation 


y@=F(t,g)+ Ly(t), q(t) + y(t), ..., On (t) + 1, v(t). 
...(4) 


Conversely, if a function y, continuous on J, fulfils (4), then the function x € C* (J) 
defined by (3) is a solution of (1) — (2). Thus the problem (1) — (2) is equivalent to 
the problem of solving integral equation (4). By substituting in the equation (4) 


eka, Tacisnia, Poze) 
= F (0,91 (t) + os 92 (t) + Piy--5Gn (14 Pua) 
we get an integral equation of the form 
YO=FGHIO, LY Os > I (O) = THO m6 
with which we shall deal. 
We make the following hypotheses used throughout this paper. 
(A,) Suppose that 


(i) there exists a continuous function g:J x Ry" > Ry = [0, o9), nondecreas- 
ing with respect to the last » variables such that 


gir, 0. 0...,,0,.0) =a0; 
itor. ihe tiecavileto (hte. Busts Furi). € Fox" 
eT ire Tin cat Tea) Ft ot Phe, Tas) | 
Set, ro— Fol, Vr Als ee | tena — Fail). 
(A,) There exists a continuous function 7: J > Ry, satisfying the inequality 


Mu (t) + A(t) < u(t) 


where 
Mu(t)=g(t, liu (t), but), ..., J,u(t) 
and 
hij sup 1 F (E,-0,.0,«, 0) | 746) 
tosESt 


(A;) In the class of functions satisfying the condition 


O<u(t)h<a(t), red 
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the function u,u(t) = 0 for t € J, is the only measurable solution of the 
equation 

u(t) = Mu(t), t€ J wr 
where Mu is defined in (6). 


In order to prove the existence of a solution of equation (5), we define the se- 
quence {)m} by the relations 


Yo (t) = 0 } (8) 

Ym4i (t) = Ty,,(t) 

fort € J and m = 0, 1, 2.,.... 

To prove the convergence of the sequence {ym} to the solution } of eqn. (5) we 

define the sequence {um} by the relation 
Uy (t) = a(t) 
Um41 (t) = Mum (t) 

fort € J and m = 0, 1, 2, .... 


gtd) 


Now we shall state our results to be proved in this paper. 


Theorem 1—Suppose that the hypotheses (A;) — (As) hold. Then there exists 
a continuous solution 9 (t), t € J of eqn. (5). The sequence {y,,} defined by (8) con- 
verges uniformly to } in J and the following estimations 


1d (t) — Yn (t) | < um (t), t E J, m= 0, 1, 2, ..., ...(10) 
and 
lp@|<a(t), tEeJ ALD 


hold. Moreover, the solution 9 of eqn. (5) is unique in the class of function satisfying 
the condition (11). 


Our next result gives conditions under which eqn. (5) has at most one solution, 
these conditions do not guarantee the existence of a solution of eqn. (5). 


Theorem 2—Let hypothesis (A,) be fulfilled. If the functionr,r(t)=0, t€ J 
is the only nonnegative, finite and measurable solutio 


r(Q@t)<Mr(t), tE J 


n of the inequality 


ati) 


then eqn. (5) has at most one solution on J. 


In order to establish our next result which de 


als with the continuous dependence 
of solutions on the right side of eqn. (5), we consid 


er the equation 


7) = KG i 2 ), hee ely (t)) -- (13) 
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where K: J X R" > Ris acontinuous function. 
Theorem 3—Assume that the hypothesis (A,) holds and 
(i) } and Zare solutions of eqns. (5) and (13) respectively; 


(ii) the sequence {v,, (t)}, t © J, defined by the relation 


MVM SIP! +124) | 

Vm+1 (t) = Mvm (t) +h (t) oP 
for t€ J, m = 0, 1, 2,..., where 

A(t) = |Tz(t) —z()| 2413) 
has a limit 2 (t) for t © J. Then 

II -ZH)|<7(), 1 EJ, ...(16) 


We note that, in the particular case where p; (t) = 1,0 <i <n, equations (1) 
and (2) reduces to 


xi*) (¢) == F(t, x (t), x’ (t), ..., x (8) Belts) 
and 
ee ig yee ya 2 =e 173. eR, ...(18) 


and consequently our results in Theorems 1-3 covers the study of equations (17)-(18). 
Here it is to be noted that the papers of Fink and Kusano’ and Trench"? are devoted 
to the study of asymptotic behaviour of solutions of (1). 


3. Proors of THEOREMS 1-3 


Before we start the proofs of Theorems |-3, we first prepare the following Lemma 
needed in our further discussion. 


Lemma -—If the condition (i) of hypothesis (A:) and hypothesis (A,)~(As) are 
satisfied, then 


0 < umy, (t) < um (t) < H(t) ..(19) 
fort € J, m = 0, 1, 2, ..., and 

Um => 0 for m — oo 
where the sign = denotes the uniform convergence in J. 


The relation (19) follows by induction. Since the sequence of continuous func- 
tions um is nonincreasing and bounded below, it is convergent to a certain measurable 
function ¢ such that 0 < ¢(t) < a(t)fort € J. By the Lebesgue theorem and the 
continuity of g it follows that the function ¢ satisfies equation (7) and by assumption 
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(As) we have ¢ (t)= 0, t € J. The uniform convergence of {u,,} in J follows from 
the Dini theorem. This completes the proof of Lemma. 


In order to prove Theorem 1, first we prove that the sequence {yn (t)}}, tE J 
satisfies the condition 


| ym (#1) | < a(t), ¢E J, m= 0,1, 2, .... ...(20) 
Obviously 
ly |= O< a(t), 6 € J. 


Furthermore, if we suppose that the inequality (2) is true for m > 0, then by the de- 
finition of y,, (t), t © J and by hypotheses (A;) and (Az), we have 


<M | ym(t)| + h(t) 
Mia (t) + h(t) 

i (t) 

fort € J. The relation (2) follows by induction. 


Ym) (1) | 


I 


In 


Next we prove that 


| Ymtq (t) — Ym (t) | <-u,, (t), t E J, m, q@ = 0, 1, 2, .... hos 
By (20), we have 
| Yq (t) — Yo (t) | = | y(t) | 


S u(t) = uw (t) 
fort € J, q = 0, 1, 2, ..., Suppose that (21) is true for m,gq > 0, then 
| Ymagts (t) — Yonr (t) | = | Tyme (t) — Tym (t) | 
SM | Yngg (t) — Ym (t) | 
S Muy, (t) = um, (t). 


Now we obtain (21) by induction. Because of Lemma Um (t) > Oin J, we have from 
(21) y,, + in J. The continuity of } follows from thé uniform convergence of the 
sequence {)} and from the continuity of allfunctions y,,. If q — co, then (21) gives 


estimation (10) and the estimation (11) is implied by (20). Itis obvious that Dis a 
solution of (5). 


+ : ; 
O prove that the solution 9 of (5) is unique, let us suppose that there exists 


another solution y of (5) such that } (¢) = y(t) and| y(t)| <a (t)t € J. By induc- 
tion we get 


A 
ly @ — yn (t)| <u, (2),1€ J me OT 2 
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F A 
and hence it follows that } (t) = y(t), t © J. This contradiction proves the unique- 
ness of } in the class of functions satisfying relation (11). This completes the proof of 
Theorem |. 


To prove Theorem 2, let us suppose that there exist two solutions } and , of equ- 
ation (5) in J, 9 (t) # y (t), t © J. Now, from hypothesis (A,) we have fort € J. 


| 3d (t) ~y(l< M | j(t) — y(t) |. ...(22) 
Putting in (22), r a |» (t)— »y (1) |, t © J, we infer from (12) that r (t) = 0 for 
t€ J,i.e. p(t) =y(t),t € J. This contradiction completes the proof of Theorem 2. 

To prove Theorem 3, let 

AN} i BC Beet 4 OR Pe TACK) 
then we have 

v(t)< | Ty (t) — Tz(t) | + | TZ) — Z(t) | 

<M\|p(t)—Z(t)1 + A(t) 


= Mv(t) + A(t). ...(24) 
From (23) and (14) we observe that | 
V(t) S|IMI+IZMO| <0), 1 € J. (25) 


Now by induction, we get 
v (1) < v, (t), t € J, m = 0,1, 2, ... ase) 
Inequality (!6) is implied by (26) as m — 09, This completes the proof of Theorem 3. 


In concluding this paper, we note that the results obtained for eqns. (1) and (2) 
in Theorems !-3 can be very easily extended for the integrodifferential equation of the 
form 


L, x(t) = F(t, Lo x(t), Ly x (t), +0) Ln-1 ¥ (O), 


f H{t, s, Lp x(s), Li x (s), «.- Ln-1 x (s)] 4s) ess 


with the given initial conditions (2), where H : 12x R° > R, F: Ix R"™ —>Rare 
continuous functions. Some results concering the existence, uniqueness and asymptotic 
behaviour of the solutions of the special versions of (27) — (2) when pi (t) = | have 
been obtained by the authors®® by using different methods. The precise formulation 
of the results similar to that given in Theorems !-3 for equations (27)—(2) are quite 
straight-forward and hence we do not discuss the details of these results. 
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FIXED POINT ITERATIONS FOR NONLINEAR HAMMERSTEIN 
EQUATION INVOLVING NONEXPANSIVE AND ACCRETIVE MAPPINGS 


C. E. CuipumMr* 


Department of Mathematics, University of Nigeria, Nsukka, Nigeria 
(Received 21 January 1988; after revision 14 June 1988) 


The solution of the nonlinear Hammerstein operator equation x + KNx = f, 
where K and WN are ‘nonexpansive’ and ‘accretive’ mappings, and K also sat- 
isfies a monotonicity condition is approximated in a Hilbert space by means 
of fixed point iteration processes. 


1. INTRODUCTION AND PRELIMINARIES 
Let ¥ be a real normed linear space. A mapping T with domain D (T) and range 
R(T) in X is called monotone if for each x, yin D (T) and some real number ¢ > 0, 
the following inequality is satisfied : 


Ix — yll < Ilx — y + 4 (Tx — Ty). = (1) 


Mappings satisfying (1) for all t > 0 are sometimes referred to as accretive’. If X isa 
Hilbert space, the accretive condition (1) reduces to 


Re < Tx — Ty,x -y>2 0 ...(2) 


for all x, yin X. The accretive operators were introduced by Browder’ and Kato"*. An 
early fundamental result in the theory of accretive operators, due to Browder, states 
that the initial value problem 


il + Tu = 0, u (0) = us ad 
dt 

is solvable if T is locally Lipschitzian and accretive on X. Browder also proved that 
if T: X > X is locally Lipschitzian and accretive then T is m-accretive, i.e., the map 
(I +T), where / denotes the identity map of X, is surjective. This result was subse- 
quently generalized by Martin’® to continuous accretive operators. If H is a Hilbert 


space, Zarantonello*® proved that the operator equation 
x+Tx=h ...(4) 


for each h € H, has a unique solution provided T is monotone and Lipschitzian. 





*This research was supported by a grant from the Third World Academy of Science (TWASRG 
86-/43), L.C.T.P., Trieste, Italy. 
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For a Banach space X we shall denote by J the normalized duality map from X 
to 2** given by 

| Ix = (f* © X*: fl? = lal? = < x, f*>} 
where X* denotes the dual space of X and < , > denotes the generalized duality pair- 
ing. It is well known that if X* is strictly convex, then J is single-valued, and if X* is 
uniformly convex, then J is uniformly continuous on bounded sets’. 

A mapping 4: H > 2# with domain D (A) in a Hilbert space H is called hemi- 
continuous at x)» € D (A), if, for any x € H such that x9 + tx € D(A) for Oct 
< ax with a, > 0, and for any sequence f, > 0 withO < ¢, < ax, we have A (xo+1,x) 


w 


w 4 ° . ° . 
—» Axo, where > denotes weak convergence. A is called hemicontinuous if it is hemi- 
continuous at every X» € D(A). It is easily seen that linear maps as well as continuous 


maps are hemicontinuous. 


In the sequel we shall be concerned with operators of the Hammerstein type, 1. ¢., 
operators of the form / + A&. These operators play a crucial role in the study of feed- 
back systems (see e. g., Dolezal’, Chapter 4) and have been studied by several 
authors!®?528, [mn Sh-Chepanovich*> the following result appears : 


Theorem Sh-Chepanovich**—Let X be a separable reflexive Banach space and let 
X* denote it dual space. Let. 


(a) N:X* —+ X bea hemicontinuous monotone map ; 
(b) K: X + X* bea linear monotone map such that for some » > O and each 
uc X*, 
<Ku,u> > p ||Kul|? mr &)) 
where < , > denotes the duality pairing. Then the operator equation 
u-+ KNu =f .-.(6) 
has a unique solution for each f € X*. 


In section 2 we examine two fixed point iteration- schemes and apply them (in 
section 3) to the iterative approximation of the solution of eqn. (6). In particular, we 
shall prove that both iteration schemes converge weakly to the solution of eqn. (6). We 
conclude with an open question. 


2. Two Fixep Point ITERATION METHODS 


In this section we describe two fixed point iteration methods given by the fol- 
lowing : 


(a) The Ishikawa Iteration Process'**4 defined as follows : For K a convex sub- 


set of a Banach ; . : : 
ch space X, and T a mapping of K into itself, the sequence {x,}"2 in K 
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is defined by 


MEK att 
Xn+1 = (1 a On) Xn s by Typ ...(8) 
Yn = (1 — 6.) Xn + B, Tx,, 2 DO ---(9) 


where {2,}°2, {B8,)°2,, satisfy0 <a, <8, <1 forall a, 


lim 8, = 0; and 2 a, B, = ©. 


(b) The Mann Iteration Process'”** which is similar to the Ishikawa iteration 
process above but with 8, = 0 and different conditions placed on ~,. More precisely, 
with X, K and x, as in part (a), the Mann iteration process is defined by 


xEK ...(10) 
X41 = (1 — C,) x, + C, Tx,,-n > 0 11) 
where (Coy is areal sequence satisfying Cp = 1;0 < C, < 1 forall n > 1, and 
: C, = cc. Thecondition 2 C, = co is, in some applications, replaced by 


2C, (1 — C,) = co. 


The iteration processes described in (a) and (b) above have been studied exten- 
sively by several authors and have been successfully employed to approximate the fixed 
points of several nonlinear mappings in Banach spaces (when these mappings are already 
known to have fixed points) and to approximate solutions of several nonlinear operator 
equations in Banach spaces*’®’$~13"161920-24. Tt is worth mentioning here that even 
though the iteration scheme (b) is similar to (a), the two schemes may exhibit different 
behaviours for different classes of nonlinear mappings**. 


3. WEAK CONVERGENCE OF THE HAMMERSTEIN OPERATOR 
EQUATION IN HILBRET SPACE 


We shall need the following results : 


Lemma®!—Let H be a Hilbert space and let {x,}°2, be a sequence in H such that 


{x,}°2, converges weakly to x* in H. Then the inequality 


lim inf ||x, — yl| > lim inf ||x, — x*|l 
nwo n—->co 


holds for all y ~ x*. 


Let K be a nonempty convex closed subset of a Hilbert space Hand let Tmap K 
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into K. T is called demiclosed at 0 in K if {x,}°2, is a sequence in K which converges 


co 
a= 


weakly to x* in K, and if {7x,}", converges strongly to zero, then 7x* = 0, 


A mapping T: H > H of a Hilbert space H into itself is called nonexpansive if 


7x — Tyll < |lx — »Il 


for each x, y in H. It is well known? that if 7: H > H is nonexpansive then (J — T) 
is demiclosed at 0. 


Theorem — Let H bea separable Hilbert space and let C be a nonempty bounded 
closed convex subset of H. Suppose 


(a) N:C + Cisa nonlinear nonexpansive monotone map; 
(b) K:C-— Cis a nonexpansive monotone map; 
such that for some » > 0 and each x € H. 


<Kx,x> > p ||Kxll*. 


Define S: C > C by Sx = f — KNx. Let {x,}°2, be a sequence defined iteratively 
by % € C, 
Xn = (1 — ay) X, + ay Syp site) 
Van CL = Ba) Xe By iSX ao oe 0) ...(13) 
where {a,} and {8,} are real sequences satisfying the following conditions : 
(i) 0 «,,8, < 1 for all n; 
(ii) lim sup 8, < 1 


iy SG poten 


Then {x,}%2, converges weakly to the unique solution of 


x + KNx =f. ..(14) 


| Proor : We shall make use of the following inequality which is valid in every 
Hilbert space, H. For each x, y, zin H, and each real number A € (0, 1), 


Ax + (1 — A) y — 2]? = A fix — aye 4 (1 — A) |ly — 2]? -a qa — a) 
lx — yl. nats) 


implies its hemicontinuity. So, the existence of 
heorem Sh. Let q denote this solution. Observe 


Observe that the nonexpansiveness of N 
a unique solution to (14) follows from T 
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that q is a fixed point of S. Furthermore, for arbitrary u,v € H, ||Su — Sv|| < |lu—vIl. 
From (12) and (13), using (15), 


Ine — gil? = [IL — oq) X, + & Syn — all? 
= (1 — a,) lx, — gil + a, IISy, — ll — % [1 — a] 
X [Xm — S¥all? SCL = aq) [Xn — Ql? + on lyn — ll? 
ae! . (16) 
We — a = I(t — 6.) x, + 8, Sx, — qll* 
= (1 — By) Ix, — al? + By Sx, — Il?—By [1 — By] Xp — Sxall? 
S (1 Bp) [Xn — I? + Bn llXa — Ql? — Ba fl — By) lle — Sxall? 
= ix, — ai — B, [1 — §,] ix, — Sx,ll?. .-(17) 
Substitution of (17) in (16) yields 
I¥n+1 — ll’ < IIx, — gil? —®, 6, [1 — 8,] lx, — Sxall*. 
Hence, ||x,.1 — gil < |lx, — q|| and {\|x, — q\|} converges. Moreover, 
nm Ball — Ba] x, — Sxall? < IIx, — Ql? — lixne — Il. 


Summation of this inequality from | to N yields 
N 
bs a, By (1 — By] Ixy — Sxl? < |lxi — ll? — Ilxw4: — II}? < 00. ...(18) 


Now, lim sup $, < l implies for AN sufficiently large, 1 — 8; > a > 0 forall j > N 
and some fixed real number a, so that the condition 2 «, 8B, = co and inequality (18) 


now yield, lim inf ||x, — Sx,|| = 0. The boundedness of {xn} implies there exists 
n—->0co oe, 


a subsequnce {x,,}7_, of {x,}2, such that {x,,},2, converges weakly to some x*€ H. 


Moreover, {x"}7°, is in C and C is weakly closed (since it is closed and convex), so it 
follows that x* E€ C. Also, 


lim IXn, — Sx,,|| = lim |Z — S) x,,|| = 0. 
ko k 


noo 


Nonexpansiveness of S implies (J — S) is demiclosed at 0, so it follows that (7/— S)x* 
= 0, i.e., x* is a fixed point of S. By uniqueness of the fixed point, x* = g. Thus 


eo 
a= 


any weak cluster point of {x,}°2, is a fixed point of S. A standard argument’”*' using 


fo) 
n=0 


the Lemma above now shows that {x,}%, has a unique weak cluster point so that 


{x,} 2, converges weakly to q, completing the proof of the Theorem. 
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Corollary—Let H, N, C, f and S be as in the above Theorem. Let {x,}"2, be a 


a= 


sequence defined iteratively by 
X EC, 
Xexic = (1: — ay) X_ te Kp Xe Bee 0 oee(19) 
where {a,} is a real sequence satisfying 
(i) 0 & «a, <1 forall a, 
(ii) af fs Target Re 


Then {x,}°2, converges weakly to the unique solution of 


x + KNx = f. 


Proor : Set B, = 0 for alln, in equations (12) and (13); and replace the con- 
dition (iii) = «, B, = co by (iii) 2 «,(1 —a,) = co in the above Theorem. Then 


the Corollary follows immediately from the Theorem. 


Comments—If K = J (the identity map of C) the Hammertsein equation x + 
KNx = f reduces to the equation 


x + Nx & f, ... (20) 


Equation (20) has been studied by several authors®”**, Dotson® showed that if N: H 
-~ His nonexpansive and monotone, the Mann iteration process converges strongly 
to the unique solution of (20). This result was generalized by the author® to operators 
with Lipschitz constant L > 1 and to operators which need not be defined on the 
whole of H. Bruck®, considered equation (20) when T: H > His a multivalued non- 
linear monotone map and proved, without any continuity assumption on N, that the 
Mann iteration process converges strongly to a solution of equation (20), if the initial 
guess is taken in a certain neighbourhood of the solution. This result has also been ex- 
tended by the author to L» spaces for Pp > 2, The methods used earlier®’”* to esta- 
blish the strong convergence of the Mann iteration process the unique solution of eqn. 
(20) seem not to be applicable to the Hammerstein equation (14) with nonexpansive, 
monotone maps N and K. This leads naturally to the following problem: 


Problem—Does any of the Mann or Ishikawa i 
to the solution of the Hammerstein equation 
monotone? 


teration process converge strongly 
(14) when K and N are nonexpansive and 
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ON WAGNER SPACES OF W,-SCALAR CURVATURE 
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The principal purpose of the present paper is to introduce and study the 
notion of a Wagner space of /’-perpendicular scalar curvature. 


1. INTRODUCTION 


Hashiguchi and Varga! examined Wagner spaces of W-scalar curvature as a 
generalization of Finsler spaces of scalar curvature and showed a result on Wagner 
spaces similar to that of Berwald spaces which was obtained independently by 
Numata’ and Varga!', Numata? gave a generalization of the above results. 


On the other hand, Izumi and Yoshida® introduced and studied the notion of a 
space of perpendicular scalar curvature (abbreviated p-scalar curvature). Later they® 
gave correct form of the equation that characterizes a Finsler space of p-scalar 
curvature. 


The main purpose of the present paper is to generalize the notion of a space of 
p-scalar curvature to a space of W-perpendicular scalar curvature (abbreviated W,- 
scalar curvature) and to prove the main result : 


Theorem— If a Finsler space is an s-Wagner space, and of W,-scalar curvature, 
then the space is conformal to a Berwald space of p-scalar curvature. 


Throughout the paper we shall use the terminology and notations of Matsumoto’s 
monograph’ with a few changes‘. For example, Uu;) means interchange of indices | a 5 
and subtraction; oj) does cyclic permutation of indices i, J, k and summation for the 
expression in the brackets behind it, and 4 := B means that A is defined by B. 


$1. WAGNER SPACES OF W,-SCALAR CURVATURE 


An n-dimensional Finsler space F, = (M,, Z) with a Finsler metric L (x, y) is 


€ exists a vector field 5; such that Fy, of 


Os; are functions of x! only'?. If s; (x) 
‘), then F, is called an s-Wagner space. 


called a Wagner space of dimension 2, if ther 


the Wagner connection WI with respect to s; 
further satisfies s; (x) = a, (s) (@% := a/ex 
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of metric tensor g;, and put hy = gi, — 1 1, (2, := 0/dy', I; : = 8, L). 


. First we derive a Wagner analogy of the known identities*"°, That is, we shall 
obtain for an s-Wagner space : 


Rij — Riku = Oc | Cur Ruy + Cry, Ri, } ...(1.1) 


C(h jk) | Rave — Chir Rj, } = (0 ... (1.2) 


r 
where XR, are components of the (v) h-torsion tensor, Rij, ate components of the h- 


curvature tensor of WI’. To establish it, we consider the expression’ (10.18') 


i i i r 
Risk = Kit Ci. R i, PAG HEY 


where 


i r 
Ki, = Yue 13 Fi, + Fr PF } 


and 5, := @, — N. 6a. For an s-Wagner space with respect to WI noticing 


St = 5, oy F' = Fi, 


r ? ¥ c 


one can easily get (A ;x) {Kia = 0, which implies 
O¢hyk) (Kutjx} = 0. ICA) 


Then on account of (1.3) we obtain 
(ijk) 1 ee One Ri, } =0 


and thus (1.2). [(1.2) is also obtained in Hashiguchi and Varga’ identity’ (11.1’) of 
WT.) 


Since WT is h-metrical, we get Rrijx = — Rinjx and so 
Kiujx + Kink = —2Cur Ry - ...(1.5) 


Now proceeding in similar way as in Rund?® (§2, 2.25), we obtain (1.1). Here, 
on account of these identities (1.1) and (1.2), the results analogous to the Lemma 1.2°® 
hold. Therefore we can apply the corresponding results in above lemma to prove 


Theorem 1.1. 
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We now consider a tangent vector ¥ = (X‘) of F, at (x, y) and the ()h-torsion 


tensor Re of WI. The quantity K = K(x, y, X) defined by 


Rio, X* X* = KL? hy, X* X* sss( 456) 
at (x, y) is called the W-sectional curvature with respect to WT. 
Definition\— A Finsler space F, is said to be of W-scalar curvature K with res- 
pect to WI’, if the W-sectional curvature K in (1.6) is a scalar field which does not 
depend on X. 


For an s-Wagner space of W-scalar curvature K, Rio = KL* hy holds 
(Hashiguchi and Varga' and Numata’). 


Now we are concerned with two independent vectors ¥ = (X‘) and Y = (Y¥*‘) in 
F,. For the plane a (p.X, p.Y) spanned by p.X‘ and p.Y', the W-perpendicular sec- 
tional curvature (abbreviated W,-sectional curvature) R:= R (x, y, 7 (p.X, p.Y)) with 
respect to WT is defined by 


Ritje (p.X") (p.Y*) (p.X/) (p.¥*) 
~ (Bh; Bik — Bae Btj) (p-X*) (p.¥*) (p.X/) (p.Y*) 


where Ry: j« is formed from the coefficients of the Wagner connection WI. 
Thus we introduce : 


(1.7) 





Definition— A Finsler space F,, (> 3) is said to be of W,-scalar curvature R with 
respect to WT, if the W,-sectional curvature R in (1.7) is a scalar field which does not 
depend on X and Y. 


Now we can prove : 


Theorem |1.1— An s-Wagner space of W,-scalar curvature is characterized by 
P.Rije = Ucey 1k Any hin + (0%, Crik + Q', Cu,)} ...(1.8) 
where 
Q,, := PR 


Proor : Similar to the proof of the corresponding result in Izumi and Yoshida‘, 
and will be omitted. 


From (1.7), it follows that the following condition 
P.Rnijx = R (hay hie — Iyy hi;) euhl ed) 
is a sufficient condition for a Wagner space to be of W,-scalar curvature. 


So, we may give the following : 
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Definition— A Wagner space F, (n > 2) characterized by (1.9) is called a 
Wagner space of Wr,-scalar curvature. : 


§2. PROOF OF THE THEOREM 


Let a Finsler space F, = (M,, L) be a Wagner space with respect to a gradient 
vector field s; (x) = 4, 5, ic. F, is an s-Wagner space). Let L* = e-*") J bea con- 
formal transformation of Finsler metrics. 


Now we consider the Finsler space F* = (M,, L*) and define 
ES = Fy, _ 5, Sk 
NM = NL — y! Se, 
c= Cr 


Since Fi of the Wagner connection are function of x‘ alone, so are the F n of the 
Cartan connection CI’ (Hashiguchi and Varga’). Hence F* is a Berwald space. We 


denote, corresponding to quantities in F,, quantities in F* by an asterisk. 


Next we assume that F, is of W,-scalar curvature, and so (1.8) holds. 
Then since 


* _ 9-2 ¥ eed 
B45 =e * 81) hy = @ ~ hy, 


Le peas Tr * pais =a2g 
Ris = Riss PeRaae = 077 P-Rirjn 


and 
by = e-%s Gage 


rik 


from (1.8) we have 
p.Rey, = Yun 18 es h* hi, + i( ar Cri + OF cs, 


and consequently, by virtue of Theorem 1.1 of Izumi and Yoshida’, the Berwald space 


F* is of p-scalar curvature R* := Re?*, which completes the proof. 


Similar to the above theorem we have the following result : 
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‘If a Finsler space is an s-Wagner space, and of Wr,-scalar curvature, then the 
space is conformal to a Berwald space of R,-scalar curvature”’. 
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Here we prove results on Goldie near-rings which are in some sense analo- 
gous to some results on Goldie theorem on Goldie rings. In some special 
types of near rings if the near rings of quotients are Goldie then the near 
rings are Goldie. An Abelian Goldie near ring K which is semiprime with 
respect to K subset and in which the non nilpotent elements are distributives, 
has a classical near ring Q of right quotient which is right Artinian and 
possesses no nonzero nilpotent right Q-subsets. 


1. INTRODUCTION 


In this paper we introduce the notion of a Goldie near-ring and prove some 
results in some sense analogous to Goldie theorems for Goldie rings. 


A countable ordered family {A,, A,, ...} of subsets of a (right) near-ring K is an 
independent family if for alln € N, 4: oe Ax) = 0, where 1 Gi < 2” and 
i 


laken. 
A (right) near-ring K is a right Goldie near-ring if 
(1) Ksatisfies the ascending chain condition (a.c.c.) for right annihilators and 
(2) Khas no infinite independent family of nonzero right K-subsets of K. 


The ring (Z, +, .) of integers and all finite near-rings are Goldie near-rings. 


We note that a left annihilator of a subset of a (right) near-ring K is always a 
left ideal of it. Therefore a left Goldie near-ring can be defined as a right near-ring K 


with the following conditions : 
(1) K satisfies the a.c.c. for left annihilators. 
(2) Khas no direct sum of an infinite numbers of left ideals. 


In our present study we shall confine ourselves to right Goldie near-rings only. 


A K-subset P of a near-ring K is ‘prime’ if for any two K-subsets A and B of K, 
ABC K=>AC KorBCK, we shall see the existence of such a K-subsets ina 


Goldie near-ring K. 
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Here we shall call a near-ring ‘prime’ if 0 is a prime K-subset, i.e., if for any two 
K-subsets A and Bof K, AB=0=>A=O0o0r8B=0. 


A right (left) K-subset A of K is nil if each element of A is nilpotent. 


A right (left) K-subset A of K is ‘nilpotent’ if there exists n € Z* such that for 
any @, ...,a, © A, a ..., @, = 0. There exist examples of such K-subsets of a near- 
ring K. (Clay’, 2, 13). 


We shall call a near-ring K ‘semiprime’ if K has no nonzero nilpotent K-subset 
of K. 


We can prove, as in case of Goldie rings (a ring with the a. c. c. on annihilator 
ideals and having no nonzero infinite direct sum of ideals) that the direct sum of two 
Goldie near-rings 1s again Goldie. We can also prove the following interesting results: 
In a semiprime Goldie near-ring the collection G of minimal prime K-subsets of K is 


finiteand (\ P=O. Another interesting result on semiprime Goldie near-ring is 
PEF : 
that if for some P € G, P = r(A) the right annihilator of a distributive K-subset A, 


then the quotient near-ring K = K/P is a prime Goldie near-ring. 


We know that: If K is a regular near-ring whose idempotents are central then 
(K, +) is Abelian*. And also there exists Abelian near-ring whose idempotents are the 
only non-nilpotent elements?. We now consider Goldie near-rings which is additively 
commutative (Abelian near-ring) and whose non nilpotent elements are distributive 
with distributively closed right essential K-susets. In this paper we establish that if 
the classical near-ring O of right quotients of an Abelian near-ring K in which non 
nilpotent elements are distributive as above and all the idempotent elements of Q com- 
mute with the regular elements of K is such that Qisa Goldie near-ring with d.c.c. on 
right Q-subset and possesses no nonzero nil right Q-subset, then K is a Goldie near-ring 
and has no nonzero nilpotent right K-subset. Next we prove that a semiprime Abelian 
Goldie near-ring in which non-nilpotent elements are distributive as above has a 


classical near-ring O of right quotients which is right Artinian and POSsesses no non- 
zero nilpotent right Q-subsets, 


2. PRELIMINARIES 


Here we assume that a near-ring K contains unity and for everya € K,a.0. =0, 


We write QO (K) to denote the 


complete near-ring of right quotients of K 
(Barua’). 


A near-ring K satisfies the (right) Ore condition with Tespect to a subset S of it, 


if given (a, r) © Kx S, there exists a common right multiple ar’ = rq’ such that 
(a’,r') © KxS. If K satisfies the (right) Ore condition with respect to S, the set of 
nonzero divisors (regular elements) of X, then the subset O={a71€ O(K) | (@,r)E Kx S} 
is a subnear-ring of Q(K). Q is the classical] near-ring of right quotients of K, 
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Let K be a subnear-ring of a near-ring L, A any subset of ZL. Then 
rp (A) = {x € L | ax = 0, forallae A} 
rx (A) = {x © K| ax = 0, foralla€ A}. 
Clearly rz (A) is a right L-subset of L and rr (A) is a right K-subset of K. 


Forany 4,BC K,A+B={a+b|a€ Abe B}. A countable ordered 
family {A,, Ao, ...... } of subsets ofa near-ring Kis an ‘independent family’ if for all 
nE N,AiN spins Ag) = 0, wherel Ci Gnand1 <k< 7x, 

i 


Let K be a near-ring and A, B are right K-subsets of K such that A C B. Then 
A is right K-essential in B if for any nonzero right K-subset C contained in B we have 
A ()} C0. A subset A of K is ‘right essential K-subset’ if it is right K-essential in K. 


A near-ring K is ‘(right) non singular’ if for any right essential K-subset A of K 
and z € K,zA = O implies z = 0, 


A ‘right annihilator ideal’ J = r (s) (=rx (s)) of K is a right annihilator K-subset 
of K such that / is a rightideal of K. 


A near-ring X is right Artinian if it satisfies the d.c.c. on right ideals. 


Since a semiprime near-ring K can not have any nonzero nilpotent right (left) K- 
subset and since in a Goldie near-ring it is possible to choose a maximal right an- 
nihilator of the type r(a) where a is a nonzero element in a right (left) K-subset A of 
K, a semiprime Goldie near-ring K has no nonzero nil right (left) K-subset. 


We now give some results (Lemmas) for use in the proofs of the main results in 
§3. The following are easy to prove. 


Lemma 2.1.1— If sy; ..., 5, € S, the set of regular elements of K, then there 
exist k,, ...k, & K,s © S such thats,’ = k, s-1,i = 1, ..., n. 
For any subset A of K we write 


QA={iqa|qa€ O,a € A} 
fin. 


AS ={Sasz la EAE s}. 


‘fin. 
If J is a right K-subset of K then by Lemma 2.1.1 we get 
IS3 = Gs) je J, s.€ S}. 
Lemma 2.1.2— JS~™ is a right Q-subset of Q. 


Lemma 2.1.3— If K is additively commutative, then so is Q. 
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Lemma 2.1.4— If K is an Abelian near-ring and A C K then 
QAQ ={ixay| x, € O,a € A} 
fin. 


is a right Q-subgroup of Q. (The result follows from the right distributivity in Q). 
Lemma 2.1.5— If TC Kand J = r(T), then JS = rg (T) and JS? 1: K=J. 


Lemma 2,2.1— Let the non nilpotent elements in the near-ring K be distributive. 
If K satisfies the d.c.c. on right K-subsets of K then every non nil right K-subset / of K 
contains a nonzero idempotent element. 


Proor: Let F = {J C /| J is a non-nil right K-subset of K}. 


Since / € F, F ~ ¢ and thus F contains a minimal element, say /;. /, being non 
nil J," is also non nil and /,* C J;. So, by minimality of J,, 1,2 = J;. Now consider 
the family of all non nil right K-subset J of K with the property that J J, 4 0 and 
J © 1,, This family is also non empty, for /, is an element in this family. Thus it 
contains a minimal element, say J;. Thus J; 7; + Oand J, isnon nil. Let u € J, be 
anon nilpotent element. SoU /; = Oandu J, C J,. By minimality of J,, ul; = # Be 
Therefore for some a € 1, u a= u, which gives u = ua" foralln € Z+. And u being 
nonzero, a@ is non nilpotent and so it is distributive. Let A (u) = {fr € J, | ur = 0}. 
A (u) is a right K-subset of K. Andu & A (u),u € J, gives A (u) C hh. So by mini- 
mality Of /:, A (u) is nil. Since uv is non nilpotent, it is distributive. Therefore u (a*—a) 
= ua” — ua = ua — ua = 0, Thus a® — a € A (u). Hence a? — ais nilpotent and 
let forn € Zt, (a2 — a)” = 0. Since a is distributive we get on expanding, 
a"=a"*! g(a), where g (x) is a polynomial in x with coefficients + 1 or— 1. Since a® 
is also distributive, we therefore get a"g(a) = g (a) a". Now a” = a"* g(a) 
= a(a" g(a)) = a (a"* (g (a))?) = a®*2 (g(a))?. Continuing this process finally we 
get, a” = a” (g (a))’. We write e = a” g(a)". Since a” € J. Andit can be seen that 
e ~ 0 and e? = e, 


Lemma 2.2.2— If K in Lemma 2.2.1 is Abelian and K has no nonzero nil right 
K-subset, then for every nonzero right K-subgroup J of K there exists an idempotent 
element e € K such that / = eK. 


Proor : By Lemma 2.2.1, / contains an idempotent element, say e!. Now 
A (e') = {x € 7 | e' x = 0} is a right K-subset of K. The family of all such sets A(e'), 
where e? is an idempotent in /, possesses a minimal element, say A (e). If A (e) #0 
then it contains an idempotent, say €,. So ee, = 0. Write e, = e + e, — e, e. Since 
Kis Abelian and e, e; are distributive we get e,? = e,. And J, being a right K-subgroup 
e, E I, Moreover, e, x =0 implies ex = 0 for ee, = e. Therefore A (e,) C A ¢), 
Again e; E€ A (e) bute, & A (e,), fore,e, =e, ~ 0. Thus A (e,.) C A (e). And 
minimality of A(e) therefore implies that 4 (e) = 0. So for any GU exer 
implies x = 0. But for anyx € l,e (x — ex) = 0. Therefore x = ~ Hence 
I= eI C eK C J which gives / = eK. . 


GOLDIE NEAR-RINGS 145 
It is easy to see the following : 


Lemma 2.3.1— Let A, B,C be right K-subsets of a near-ring K such that 
A Cc BC CC Kand Ais right K-essential in B; B is right K-essential in C. Then A 
is right X-essential in C. 


Next we prove: 


Lemma 2.3.2— Let N be a right K-subset of K and M be such a right K-subset 
of K that M is right K-essential in N. 


If a € N, a ¥ 0, then there is a right essential K-subset L of K such that aL0 
andal C M. 


Proor: Write L = {k € K|ak € M}. Clearly L is a right K-subset of K and 
aLCM. Since N is a right K-subset of K,aK C N,aK 4 0and aKQ M40 
(Since 1 € K and M is right K-essential in M). Therefore there is some k € K such 
that ak € M, ak #4 0. SoaL + 0. Now let J be a nonzero right K-subset of K. If 
al = 0, then J C LZ and hence] ™ L 40. And if aJ + 0, aJ C Nand Mis right K- 
essential in N give that al (\ M+ 0. So for some x € J, ax 40, ax € M. This 
implies thatx € 7 L. Sincex 40,77 LZ 40. Thus L is a right essential K- 
subset of K. 


Taking N = Kand M = 4A, a right essential K-subset of K we get the 


_ Corollary 2.3.3— For anya € A,a* A= {x€ K |ax € A}isa right essen- 
tial K-subset of K. 


Lemma 2.3.4— Let A, B be right annihilator K-subsets of K with A C B and A 
be right K-essential in B. If K is right nonsingular then A = B. 


Proor: Letb € B,b #0. Since A C B and A is right K-essential in B, by 
Lemma 2.3.2 there exists a right essential K-subset L of K such that DL C A, bL + 0. 
Thus /(A) bL = 0 and K being right nonsingular we get /(A) b = 0 which gives 
b € r(I(A)). Since A is of the type r (S) for some S K,r (/(A)) = A. Hence b € A 
whence B C A. Thus A = B. 


Lemma 2.3.5— If A and B are two right K-subsets of a right nonsingular near- 
ring K such that A C Band A is right K-essential in B, then for any x € K the subset 
xA is right K-essential in xB. 


Proor: Let C C xB, C ~(0) be a right K-subset of K and letc = xb be a 
nonzero element of C. Then for b € B there is one right essential K-subset L of K 
such that bL ~ 0, bL C A (Lemma 2.3.2). Since K is right nonsingular, we therefore 
get, xbL # 0, for otherwise we shall get xb = 0. Nowxd 1) CD xA 1) xbK 3D 
xbL ~ 0, which implies that xA is right K-essential in xB. 


Lemma 2.4.1— Let K be a Goldie near-ring whose non-nilpotent elements are 
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distributive. If x € Kis such thatr(x) = 0, then xK is a right essential K-subset 
of K. 


Proor : Since r (x) = 0, x is non ni)potent and therefore it is distributive. Let 
M be a right K-subset of K such that M ™ xK = 0. Now fora fix s € Z* and for 
t<s, let 


eaEG(xM) Ax M, OO = la =O 1,5 8): 
n#i 


Then 
a= 2 x"m, = x' m, mm, E M. 
n#Ft 
1.6.5 
my = x'm, — x'm,—....—xm 
= x (x!"! m,., —x*-) m,-1 ... — m,) (since x is distributive). 


Thus, my) € M () xK which gives m, = 0. It follows that 
x) m-, —x*") m_, — ... m, = 0, for r(x) = 0 
Similarly we get m, = m, = ... = m = 0. Therefore ( = x" M)Q x! M = Ofor 
a5 


alls € Z*+ andi < s. Thus the family {M, xM, x? M, ...... } is an independent family. 
K being Goldie, there exists u € Zt+ such that x“*! M = 0. So for anyme M, 
x“*} m = 0 which gives m = 0 (Since r (x) = 0). Therefore M = 0. Thus MAxK<0 
implies M = 0. Hence xX is a right essential K-subset of K. 


Lemma 2.4,2— The right singular K-subgroup of a Goldie near-ring K is nil- 
potent. (Follows from Proposition 4.3.6 in Barua’). Therefore 


Lemma 2.4.3— A semiprime Goldie near-ring is right non singular. 


Lemma 2.4.4— If K is a semiprime Goldie near-ring where non-nilpotent 
elements are distributive then every distributively closed right essential K-subset of 
K contains a regular element. 


Proor : Let J be any distributively closed right essential K-subset of K. We first 
show that / contains an element a such that r (a) = 0, 


Since K is semiprime and / + 0, it is not nil. Let us choose a non nilpotent 
element a, € / with r (a) as large as possible (it is possible for K is Goldie). Suppose 
r (a,) #0. Then r (a,) Q 7 # 0. As above choose a non nilpotent element a€Er 
(a,) M I with r (a,) as large as possible. And Suppose r (a, + a.) + 0, Now r(a,+a,) 
(™ 1 0 and it is not nil. Choose a non nilpotent element a3 € r (a, + a.) \ J with 
r (a3) as large as possible. We see that r (4) Or(@) =r(at a2). Now we note that 


(a) 1 r (a,j C r (a, + a). First we show aK a,K=0, Letz =a,x=a,y 
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belong to the intersection. Since a, : a; are non nilpotent, by maximality of r (a,), 
r (a2) we get r (a, ) = r(a,) and r(a? ) = r(a,). Therefore a. x =a,a, y = 0, 


fora, € r(a,). Thenx Er (a) ), ie. x © r(a,). Thus z = 0, Hence a, KQa K=0. 


Now if x € r (a, + a,), then a, x = a, (—x), for a, is distributive. Since a,K(\\a,K=0, 
We get x € r(a,) 1 r (az). Hence r (a, + az) = r (a,) Or (a,). As above we see 
that ia, K, a, K a, K} is an independent family and r (a, + a, + a3) = r (a,) N r(a,) 
™ r (a3). We continue this process. Because of Goldie condition we get a non left zero 
divisor C = a, + ...... + a, i.e. r(c) = 0. Weclaim that /(c) = 0. Let xc = 0 for 
some x € K. Then xcK = 0. By Lemma 2.4.1, cK is a right essential K-subset of K. 
And K being semiprime, it is right non singular. So xcK = 0 implies x = 0. Hence 
I (c) = 0 Thus ¢c is a regular element of K. 


Lemma 2.4.5— Let K be a semiprime Abelian Goldie near-ring such that its 
non-nilpotent elements are distributive. Then K-satisfies the d.c.c. for right annihilator 
ideals. 


Proor : Let A, B be two right annihilator ideals. Then A = r (T) and B= r(S), 
S,T C K. If A C Band 4 is not right K-essential in B, then by Lemma 2.3.4. there 
exists a right K-subset P (CB) such that A () P = 0. K being semiprime Goldie, P is 
not nil. Let p be anon nilpotent element in P and so it is distributive and pK is a 
right ideal of K, for Kis Abelian. Now AM pKC A(1)} P=Oie. AQ pK = 0. If 
possible let B > A D Cbe a strictly deceending chain of right annihilator ideals 
where B = r (T), A = r(S), C = r(U). As above there is a right K subset 0 C A 
such that C () Q = 0. And as in case of P, there is a distributive element g € Q such 
C1 qK = 0. We claim that {pK, qK, C} is an independent family. For 4 ) pK = 0, 
Cf qK =0, pK C BC, qK CA. Hence pK (qK + C)C pKN(A+C) 
C pK) A=0. If qki = pk, +c € qK 1 (pK+ C), then pk, = gk; + (—c) € (qK 
+ C)MpK =0. Thus gkj = c € C 1) gK = 0. Similarly C M (pK + gK) = 0. So an 
infinite strictly descending chain of right annihilator ideals gives an infinite independent 
family of right K-subsets of K (for K is Abelian), which contradicts the Goldie charac- 
tor of K. Hence K satisfies the d.c.c. for right annihilator ideals. 


3. MAIN RESULTS 
We now give the main results. 
Theorem 3.1— If the near-ring Q is Goldie then so is K. 
Proor: Let J: C Jo C ...., where J;) = r(7i), T7 GC K be a strictly ascending 


chain of right annihilator K-subsets of K. Then by Lemma 2.1.5. 


J, SC Jo SC ...... is an ascending chain of right annihilator Q-subsets of 
Q. The near ring Q being Goldie, Jm S~' = J,4,; S? = ......, for some m € Z*. And 
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by Lemma 2.1.5 we therefore get Im = Jmy) = «..-+- Therefore K cannot have any 
infinite strictly ascending chain of right annihilator K-subsets. Next if ‘i, Ui See Me iheet 
_an independent family of right K-subsets of K, we claim that the family {J, S~*, ..., 
J, S~} of right Q-subsets of Q is independent. If for some m,]<m St, JmS™ (1) 
(<= J, S-) # 0, then we get a nonzero element /,, s = fles Aa in the intersec- 
n>=-m 


tion. By Lemma 2.1.1, we get k,, .., ks © K,s € S such thats’ = kis?,1<ict 
And because of the right distributivity in Q, 


In Kot = ZB jks =f SS 7, Ky) s—* which. gives 
sm n-=m 


n 
jn ke = 2 J, k,. Thus.J,.-\) Geo al.) = 0. 
nm nm 
This contradicts that the family {/,, ..., J:} is independent. Therefore K can not have 
an infinite independent family of right K-subsets. 


Thus K is Goldie. 


Theorem 3.2— Let K be an Abelian near-ring in which non nilpotent elements 
are distributive and the regular elements of K commute with the idempotent elements 
of its classical near-ring Q of right quotients which satisfies the d.c.c. for right Q- 
subsets and which possesses no nonzero nil right O-subset. Then K is a Goldie near- 
ring with no nonzero nilpotent right K-subset. 


Proor : By Theorem 3.1 K is Goldie. 


Now suppose N is a nilpotent right K-subset of K such that N? = 0. By Lemma 
2.1.4. ONO is a right Q-subgroup of Q. So by Lemma 2.2.2, ONO = eO for some 
nonzero idempotent e € Q. Thereforee = X11 Yi, X1, Ys © O, ni E N and each 

fin, 


wok s',k€ Ky ES. Using Lemma 2.1.1. We get um € K, s © S such that 
each s; = uy, s,i = 1,......, t. Therefore e = ( S ximkiuj)s. And this gives 
n. 
es = 2 xin kj us © QN (since each n; ki us € N).SoesN C ON?, i.e., esN = 0, 
in 


Since e commutes with s, seN = 0 which gives eN = 0, And N C ONO = €O gives 
that form € N,n = eq,q € QO. Hencen = en € eN. Thus N C eN, i.e., N = O, 


Theorem 3.3— A semiprime Abelian Goldie near-ring K in which non_ nilpotent 
elements are distributive with distributively closed right essential] K-subsets, has a 
classical near-ring Q of right quotients which has no nilpotent right Q-subsets. 


Proor : Choose a,b € K,a is regular in K. K being Goldie, by Lemma 2.4.1 
a Kis right essential in K. Then, by Corollary 2.3.3., the set A ={kE K|bke aK} 
is right essential in K Therefore by Lemma 2.4.4, A contins a regular element, say a’. 
Thus bat = ak!, for some k1 € K. Thus the right Ore condition with respect to the 
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set S of regular elements of K is satisfied in K. So by Lemma 5.4.4, (Barua!) Khas a 
classical] uear-ring of right quotients of K, say Q. 


Next let J be a right Q-subset of Q such that J? = 0. NowJ () K is a right K- 
subset of K. Because of Lemma 2.1.1, (J ( K) Q = {xc+ |x € J K,c € S}. Now 
for any x € J,x =ks",k € Kis € S. Soxs mk Ee JK. Thus x = (xs) SOE 
(J K) Q. Conversely if ys? € (J 1 K)O,y EC JN K,s € S, then ys! € J. Hence 
J =(J 1 K) Q. Again J? = 0 gives (J ) K*)(C J?) = 0. Therefore J 1 K isa nil- 
potent right K-subset of K and K is semiprime. Hence J () K = 0. Thus it follows 
from what we have showed above that J = 0. 


Theorem 3.4— A semiprime Abelian Goldie near-ring K in which non-nilpotent 
elements are distributive has a classical near-ring Q of right quotients which is right 
Artinian. 


Proor : First we show that if A, B are two right ideals of Q with B C A and 
BQ Kis right K-essential in A () K then A = B. Letx € AQ K. By Lemma 2.3.2, 
there is a right essential K-subset ZL of K such that xL € BQ K. Andby Lemma 
244, F contains a regular element, say c € L. Thusxc € BQ Kandx = (xc)c'E 
((B 1 K) Q =) B. Therefore A M K € 8B. Hence A = (A —) K) OQ C BQ C B. Thus 
A = B. 

So if B C A, then BQ Kis not right K-essential in A ( K which implies that 
there is a nonzero right K-subset X¥ of K contained in A M K such that ¥(BK)=0. 
Since X cannot be nil, it contains a distributive element, say x. Again xK(\(BNK)=0 
and xK is a right ideal of K. Further we have xK C X CA K. If ADB D Sa hy 25: 
where C, D are right ideals of Q, then in like manner we get right K-subsets, Y, Z of 
K contained in BOK and CK respectively such thatyK CYC BN K, 
ZERO ZOCCHK,yKO(CN K)=0 andzKN (DN K) = 0. We show that 
{xK, yK, zK} is an independent family. Let xk, = yk, + zk, © xK 1 (yK + 2K). 
Then we get yk. + zk; € xK 11 (BO K) (since BN K is a right ideal of K). Thus 
xk, = 0, ie. xK M (yK + 2K) = O. Similarly yK (\ (xK + zK) = 0 = 2K 
(xK + yK). Therefore the family {xK, xK, zK} is independent. Thus an infinite 
strictly descending chain .of right ideals of Q gives an infinite independent family of 
right K-subsets of K which is in contradiction with the Goldie character of K. Hence 


Q must be right Artinian. 
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MATRIX TRANSFORMATIONS OF ORTHONORMAL SERIES 
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(Received 6 March 1988; after revision 4 July 1988) 


Let D denote the set of Hausdorff matrices for which max, h,, =O(n71/°), E, 

k-1 

the set of lower triangular matrices with row sums one satisfying 2: 
j= 


@ny = O (k/n), uniformly in k. In this Paper we establish number of theorems 
involving the summability of orthonormal series by matrices in either class 
D or E. These results significantly extend some of the corresponding theorems 
established by Meder4 and Patel5-7 for the Euler matrix of order li 


A Hausdorff matrix H = (A,x) is a lower triangular matrix whose nonzero entries 


are of the form /n, = ( A ) A"-* 4, where : ) is the ordinary binomial coefficient, 


{u,} is a real or complex sequence, an A is the forward difference operator defined by 
A® we = pe, Aue = py — Pkziy APF? ny = A (An ax). Examples of Hausdorff matrices 
are the Cesaro matrices of order “, obtained by Setting 1, = 1/ ear ), and the 
Euler matrices, (E, p), obtained by setting wy, = (1 + Py iol 


Let » (x) be a Positive, bounded, monotone increasing function defined Over an 
interval [a, b], whose derivative is no 


Se @) = Ea $i (x) 


where the a; are real and satisfy 


> 2 < 00, -..(2) 


n=0 


Let {/ = x = {x4}: Ze | xp | <ce}; 


Then / is called the set of absolute conver- 
gent sequences, where the convergence of 5 


| x, | means the absolute convergence of 
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the series = x, rather than that of the sequence {x,}. B&(l) will denote the set of 
bounded linear operators on /; i. e., a matrix 4 = (a,x) © B(l) if, for each x € |, 


AX = {2 Gyx Xk} € 1. A Hausdroff matrix H € B (1) iff fi | dB (t) | < co, where 
0 


1 
B(t) € BV [0, l] and p, = § t" d8 (t) (see Rhoades*). Let D denote the set of Haus- 
0 


dorff matrices for which max, (A,x) = O (n-/*), 


For H a Hausdorff matrix define +, (x) = Zio Anak Sk (x). 


Theorem 1—The series 


Bf meee (2) — toa (a) du @) 3) 


is convergent if H € B(I) (\ Dand 


: S Vn a. < co, ...(4) 


Proor : It has been shown’ that, if a series Sa, has partial sums t,, then 
k 


1 S } : 
t, (X) — T.-1 (x) = a Tn aj. 
i=0 


Using the orthonormaly of {¢, (x)}, 


6 
f n [t, (x) — Th-1 (x)? dp. (x) 


b n n 
am | ae oS h,; ja; 4; (x)) (> hyx kay bx (x) du (x) 


n 
a thoes, k=O 


l . 2 2 
soos > | ee | aa, 
n nk k 


k=O 


I 


< * nl max. | Ay. | Sy | Hak | Vka} 


k=0 


ll 


O01) S| lnk | Vat 5) 


k=0 


152 B. E. RHOADES 


Since H € B (I), > >» | Ank | / ka; is finite, and Theorem | is proved. 
1 k=l 


There are many Hausdorff matrices that belong to D. We shall list here two 
examples. 
For the gamma methods, p, = a/(n + a), a real, and 


_ Tat+prk+a 
a Pk + 1H + a +1) 


har k+1 es 


h 


— 





ea ka: (21 for a ern 
Therefore 
max, hes = | hn, es a < l 
ae bo oe be 


( T@rmt)) 





Fon ie pcs ee | 
a ak Teh Perr 
{ n+a ’ i 


The mass functions for the gamma methods are @ (t) = f*. Therefore each opera- 
tor method belongs to D fora > 1/2, and belongs to B(/) fora > 1. 


For the Euler methods, using either Theorem 138 of Hardy? or Lemma 1 of 
Ziza', max, h,, = O(n-'/*). That each Euler method is in B (1) follows from 
Theorem | of Rhoades®. 


Meder‘ has shown that (3) implies (4) for (E, 1). On the other hand, Ziza™ has 
shown that all Euler matrices are equivalent a. e. for every orthonormal series with 
coeflicients satisfying (2). Therefore a reasonable conjecture is that Meder’s result can 
be extended to (E, p) for p > 0. 


Corollary 1\—Let (E, p) be an Euler matrix of order p > 0. Then (3) converges 
a. e. iff (4) converges. 


That (4) implies (3) follows from Theorem 1, since every Euler matrix belongs to 


To prove the converse, note that one can replace the convergence of (3) by that 
of (5). For N sufficiently large, and for ni(p+1)—VWn> Ini/(pt+lI)t+ayn<n 
for alln > N, condition (5) implies that 


I 
‘ SS k? h* a; — oO. 
n 


n=N n 
pt+i1 ~V"Sks pth we 
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But, from Theorem 138 of Hardy*, h,« > c/\/k over the range of the inner sum, 
where c is some positive constant. Therefore we have 


io 2) 
l ° 
$i: Mt <= 
a ey LIS 
Pel ~V8FS 54 tvs 
which in turn implies that 


co 
S S a < 9, 


laa n n 
es ae VnSk a +Vk 





The above sum can be rewritten in the form 


> N (k) aj 


where N (k) denotes the number of integers n > N satisfying n/(p +1) -—-Vn< k 
<n] (p + 1) + Vk. It can be shown that N (k) is asymptotic to Vk, and the result 
is proved. 


Let E denote the set of regular lower triangular matrices, with row sums one, 
satisfying 


"Say = O (kjn) 


j=0 
uniformly in k. We shall show that E is a rather large class. 

First, observe that, since each A in E is regular, each a,x is bounded. Therefore, 
if » is any real number satisfying 0 < 4 < 1, then, automatically, 


= a, = O (k[n) 


k > yn. Thus, it is enough to verify the equality for k < yn. 
Let (N, p) be any regular Norlund method satisfying np, = Ott Pa) ne 


n 
entries of a regular Norlund matrix are ¢nk = Pu-k/Pns Pn = Re Dk. 


k-1 
2 cpt — PH F,. 


j=o 
Claim: P,-x/P, = 1 + O (kjn). 
Pi-il fs =1-—- PalPrn =I1+ O (1/n). 
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Then log | P,-,/P, | = O(1/n) and log | P,-«/P, | = O(k/n). Therefore | P,-«/P, | 
= exp (O (k/n)) = 1 + O (kjn) 


For a regular weighted mean method (WN, p), the entries are cnx = Px/P,, and 


k-1 k=1 
Z Cy = = pjlPy = Pr-1]P, = O(kjn). 
j™=0 


j=0 


uniformly for 1 < k qn for P, ~ n*,a > 1. 


For any regular Hausdorff matrix H, {»,} has the representation 
1 
Hn = J Be) 
where 8B € BV [0, 1]. Suppose that H also satisfies 
t 
J14B(t)|= ot), ast>0 +. .. (6) 


1 
Since J | dB (t) | < ©», it follows that, if (5) is true for small values of t, then it con- 


tinues to remain true, possibly with a different constant, uniformly over O<1t< 1. 
It is sufficient to show that 


S hy; = S | & Je — 1'as (t) 
| [S ( : axe = om |e (t) = O (kjn) ...(7) 


uniformly for 1 qk & n/4. 


Since the expression in brackets in (7) is bounded above by one; it follows from 
(6) that 


2k/n — k—1 
[ oy (3 )va- 4] dp (t) = O (kjn). 
0 j=0 
It remains to estimate the contribution of the integral in (7) over the interval [2k/ n, 1] 


Since [2j/n,1] includes [2k/n,1] for 0 < J Sk — 1, itis sufficient to show that 
uniformly for 0 <j < n/4, 
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i) (1 = 2-4 | dB (t) | = O(IIn. oe) 


2jin 


t 
With X (t) = J | dB (u) |, X (1) = O(1), and the left hand side of (8) can be written 
in the form 


1 


-(" | me) — x ily) WH — D9 at. 


2i/n 


Since d (t’ (1 — t)"-’)/dt < 0 in the interval of integration, the above expression 
is dominated by a constant times 


1 


San i < (t* (1 — t)"-4) dt 


2jin 


( ) ate Plestiseoie ee ( 1 (1 — 1)" de} 


24[n 


= 1,, + i, say. 
1 
hos (ee) | (1 — t)* dt = I]/(n + 1). 
0 


For j = 0,4, = 1. Forj 2 1, 
hh = (7) Gamer (1 — 2m 
MNase: Ue 2p 
oa ee 


yntiie Qitl y2y? (n = 4 Bvt 
~ FPR (a — jr * meee 


Pei in nig gene ocean 
a eee Te 

The Hausdorff matrices with 8 (t) = 1%, 0 << «<1 donot satisfy (6) and do 
not belong to E. Therefore condition (6) cannot be weakened. 


Let o, (x) denote the nth term of the (C, 1) transform of. G)s4. ¢., 


o, (x) = (2 + 1)-7 = Sk (x). 
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The next two theorems require the following lemmas. 


Lemma 1—If the orthonormal series (1) satis fies (2), then, for any A € E with 


t, (x) = Zz Onk Sk (X) 


= [o, (x) — (Opin (9) 


converges a.e. 


Proor: o, (x) — t, (x) = 5 es — Ank | Sx (x) 


2 > i =n | S140 








; > ern > lar ~ | 
a > aj; (x) [[(n—j+1)(n+1) — > ax] 
-> asfi- — LS aw] 
oa k=j 
@1 2 a); (*) | > nk — j[(n + )] .-(10) 
Since A € E, 
co 6b 
aS | ([on (x) — ty (x)}2/n) du (x) = O(1) Sa a > a? (kjn)® 
= 00) S ee =. 
2 es o>, a 


Using t i i 
ing the theorem of Levi (see, e.g. Alexits?, p. 11) (9) converges a. e. 
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Lemma 2—If the orthonormal series (1) satisfies (2) and is summable (A) a.e. to 
s(x), for A € E, then 


ie 
lim a4 S [Sk (x) — S (x)? = 0. 


PROOF : s, (x) — o, (x) = ean > k ax ox (x). 


Therefore 


co 6 co 


> | ~ [s, (x) — ¢, (x)}? dp (x) = aS ke a? 
mr hoes = pa 
=> ae roo S at 


k=} n=k 


From the theorem of Levi, 
co 
= [s, (x) — o, (x)}/n 
n=1 


converges a.e. From Kronecker’s Theorem 


4 BS [sz (x) — o% (x)]? ~ 0 asn > co, 


k=1 


Using Minkowski’s inequality, 


> SS [x (x) — spe <{t Ss [x () — % GP pe 
k=l k=1 
+{ Ps (a (2) — 1 Gp be 
+15 S [i @) - s QP be. 


Since t;, (x) > s (x), the third series on the right converges to zero. Applying 
Kronecker’s Theorem to (9) yields 
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y) [o, (x) — t (x)]? > 0 asn > oo. 


k=l 


Theorem 2—If the orthonormal series (1) satisfies (2) and is summable (A) a.e. to 
s (x), A € E, then it is summable (C, 1) a. e. to the same sum. 


PROOF : 


[@) — s@OF <{ py, 1e@ — 51 } 


k=a0 


k=0 


M\ 


and the result follows by Lemma 2. 


Theorem 3—If the orthonormal series (1) satisfies 


co 


DS a* (log log n)? < co 


Neo2 
and 4 € E, then lim, t,, (x) exists a.e. 


Proor : From the proof of Lemma 1, 


= j [o.” (x) — 1,” (x)}? dy (x) 


I 


co on k-1 
2 n k \ 
BIS | > 3, — nee 
n=1 k= j=) 


. 


n 
co 


O (1) > 4-" Ss a’ e=ow > ay 


n=l k=] 


n=} 


by the argument of Alexits! (p. 145). 


The remainder of the proof is the same as that in Meder‘, 


Theorem 4—Let A bea nonnegative Hausdorff matrix, AE DQ) Elf the ortho- 
normal series (1) is (4) summable a.e. to s 


(x) and if (4) is Satisfied, then it is strongly 
A-summable, with index 2, a.e. to s (x). 
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Proor: © One [Sz (x) — s(x) <2 ¥ agelse (x) — te OOP 
=0 2.0 


+25 on lte (xs) — sd} =S, + S, say. 
k=0 


By hypothesis S, — 0. 
S,=2 > Onk [Sk (x) — te (XP = ons [sx (x) — te COF 
k=0 


k=0 


k 
Sx (x) — te (x) = Sx (x) — Ree ak; S; (x) 
k k j 
= = a;¢;(x)— 2% ax B a d(x) 
40 3=0 1=0 


k k k 
= ZX aej(x)- LS agi(x) & a; 
j=0 i-0 jal 


k k 
2 ay gi (x) [1 — eo ax] 


k i-] 
= D a,gj(x) = apy. 
{=0 j=0 





co [oe] n fe | 
Sa [s () — 4 PAO) = > Tae a; [> aus | 
=O o> 2S aj (i?[n*) 


co we) 
= O(l) > th Sy n~5/? 
{=1 n=4 


ioe) 
=O(1) = Via. 
f=) 


By Levi’s Theorem 


ioe) 


> si [s, Oe (PF 


n=l 
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converges a.e. By Kronecker’s Theorem, 
) . 
© [x () — OP = O (vn) 


and S, > 0. 


Theorem 5—Let the orthonormal series (1) satisfy 


co 
>> Am < co tg 9 


m=9@ 


where 


m+1 
2 
: | 2 
imo 4y 
Then 
co 


| Sn) = ty (X) | [1 < 20 


for each A in E 


PROOF : 
n k-y 
Sp (x) = ie (x) = ak dk (x) Be je 
Therefore 
ya 5 +[3 
— | 1 Sn (x) — th (x) | dp (x) = O(1) ae [ kay] 2 
n=l % a x Sy k 


Without loss of generality we may assume a, = 0. It is well known that 


n [logn] 
> ea < Sanna, 
k=z r=a0 
Thus 
co b n Crogan] 
3 Fal 156) =) | de @) < 0”) Sn > 24, 
2 a n=2 r=9 


(equation continued on P. 161) 


MATRIX TRANSFORMATIONS 161 


Cie s 2 oat as 
r=2 log[n]?r 


foe) 
OT Ae 
rem? 
Theorem 6—If {Am} satisfies (11), then, for each H in D. 
= | te (X) — Tey (x) | = 0 (on). 


PROOF : We may assume ao = a; = 0. Using (5), 


p ; | | tm (x) — toa (x) | dus (x) 





Jn 


n=l 


6 


<> A[fis@-n ere] 
Sen {Srn ap 


neal o=2 


co 


= 00 S02 {S ra 


n=1 j=2 


rere) 
O(1) = 2i+1 A; >» na 
j=2 logiut+1]7 7 


co 
= 0 (I) % A) <0. 
=2 


The proof is completed by using the theorems of Levi and Kronecker. 


We now prove two theorems dealing with lacunary series. 
Theorem 7—If {a,} satisfies (2) and {n,} is an increasing sequence of indices 
satisfying 
1<q < Mryilm, k = 0, 1, 2, Ale) 
then, for each A € E, 


E fo) — te OP a 
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converges a.e. in [a, 5]. 
PROOF : S, (x) — t, (x) = 2 ax bk (Xx) a Ggf. 
=—( = 


n 


1 |, (x) — tn (x) | du (X) |: =O (7 aS: a, (k/n)* 





| hss! 
since A € E. 
Therefore 
oo «6b ps ny 
| | Sup — tnx (X) | *du(x) = O (1) ys ja 
k=l a Sane 


ll 
So 
M 
oo 
5, 
Ms 
a 





ll 
eo) 
ae 
— 
~~ 
nw 
=) 
“™ 1 
vad Wa 
iL) 
Ri 
3 


t Aimy | m=0 
co 
= O(1) a; < 00 
j= 


Now apply Levi’s Theorem. 


| Theorem 8—If {Gn} satisfies (4) and {n,} is an increasing sequence of indices 
satisfying 


1<q4  Megilm <r, k = 0, 1, 2, ... ...(14) 


then (1) is (E, p) summable ae. in [a, b) iff {s,, (x)} converges a.e. in [a, b]. 


PROOF : Suppose (1) is summable (E,p). Then (14) implies (12) and (13) is 
satisfied. Since {r, (x)} converges, {s,, (x)} converges a.e. 


Suppose {s,, (x)} converges a.e. for {n;} satisfying (12). Let 5 Satisfy ms - 
+= Nk+}. 


SO), Ob 3 lee GemananE 


(equation continued on Pp. 163) 
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n 


"k+1 k+l 
oe peteeccR (re tv uesgee os (x))h oer Al nc nia C15) 
at 


n ann, +1 
But 


n 
k+1 


| 
IJn< a (M1e4; — Me) = Ney |me — 1 <r — I. 


otriad 2 


From Theorem |, (5) implies (2). By Levi’s Theorem 
co 
aa n [t, (x) -— Tn-; (x)]? 


converges a.e. Therefore the right-hand side of (15) tends to zero a.e. From Theorem 7 
{t,, (*)} converges a.e. Therefore {t, (x)} converges from (15). 


The author is indebeted to Brian Kuttner for the idea of extending the results of 
Meder? and Patel®~? to the classes D and E, and for the proof that those Hausdorff 
matrices satisfying (6) belong to E. 


REMARKS 


1. Theorems |—3 are generalizations of Theorems !,2, and 4, respectively, of Meder‘. 
2. Lemmas | and 2 are generalizations of Lemmas | and 2 of Meder‘. 

3. Theorem 4 generalizes the Theorem in Patel’. 

4. Theorem 5 generalizes the corresponding result in Patel®. 

5. Theorem 6 generalizes corresponding Theorem of Patel’: 

6. Theorems 7 and 8 are generalizations of Theorems | and 2 of Ziza?® 


7. In Sharma® the result of Patel® is extended to double Euler summability (E, 1, 1). 
However, the results in that paper are incorrect since the author assumes that 


m n m m #1 /-1 

>» x— 2 = > ae 

j=0 k=0 77 k=l j=0 k=0 
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STRESSES IN PRE-STRESSED DRY SANDY SOIL DUE TO NORMAL 
MOVING LOAD LEADING TO INSTABILITY AND FRACTURE 


S. Dey AND M. CHAKRABORTY 
Department of Applied Sciences, Indian School of Mines, Dhanbad 826004 


(Received 25 August 1987; after revision 4 July 1988) 


The paper gives a complete study of stresses and displacements induced ina 
pre-stressed half-space made of dry sandy soil due to concentrated line load 
moving at a constant speed along the surface. The subsonic, transonic and 
supersonic cases have been considered. It is seen that sandy parameters 
and the pre-stress parameters play important roles in the development of the 
Stresses and displacements in the medium. The stress developed at a certain 
depth have been calculated numerically for the subsonic case due to increas- 
ing velocity of the moving load for different values of pre-stress parameters. 
The velocities of the moving load creating instability in the medium leading 
to fracture have been calculated for different values of pre-stressing in sandy 
soil and elastic half-space both and it is observed that sandy soil is less stable 
to moving loads than the elastic one. Further it is inferred that the pre- 
stressed medium gets fractured at a less velocity of the moving load in com- 
parison to the pre-stress free case. 


INTRODUCTION 


To analyse the stresses developed in a body due to a moving source causing frac- 
ture is an interesting problem of mechanics having its application towards the stability 
of a medium. Sneddon'” has developed an analysis which gives us the displacement 
on the surface of a semi-infinite medium for different kinds of source of disturbance 
applied on the surface. Of these, the particular kind of source which is acting parallel 
to the surface and moving with a certain uniform velocity is of special interest to seis- 
mologists. In the above paper it has been shown that when the velocity of moving load 
exceeds the velocity of shear waves inthe medium, displacement becomes infinitely 
large along two lines. The steady state solution of the problem of moving normal load 
over an elastic half space were given by Cole and Huth* and Craggs*, who derived a 
relatively simple closed form solution, exhibiting a resonance effect at a critical load 
velocity, which in this case equals to the velocity of Rayleigh waves. The problem 
considered by Cole and Huth® has been discussed previously by Sneddon’ by a some- 
what different method. However, Sneddon! treated only the subsonic case. Ghosh® 
has explained the principle behind the phenomenon of propagation of cracking across 
the length on the simplified assumption that the crust is lying over a medium (rigid 
foundation) with shearless contact and a normal point source is moving with a certain 
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velocity. Stresses developed in a transversely isotropic elastic media (under different 
conditions) due to normal moving load over a rough surface have been discussed by 
Mukhopadhyay*, Mukherjee? and Dey, ef a/.*. Freund? discussed wave motion 
in an elastic half-space subjected to non-uniformly moving surface load. 


In this paper attempt has been made to study the stresses and the displacements 
developed leading to fracture of more earthy material, say, dry sandy material under 
initial stresses. The crust of the earth is not exactly elastic but may be estimated as 
sandy material whose definition was given by Weiskopf?®. Further, the normal initial 
stresses developed in the earth due to many physical causes such as variation of gravity, 
temperature, slow process of creep etc. deserve its consideration in development of 
Stresses and displacements due to a moving load on the surface. 


The relation E/» = 2(1 + s) for isotropic elastic solids does not hold good for 
real earthy materials viz. sand, soil etc. Weiskopf'® investigated that due to slipping 
of granules on each other the resistance of shear is much less than that ina solid and 
the resultant shearing deflection is much greater. For these materials 


E/up > 2(1 + a). 


So the relation E/ = 27 (1 + a), where 4 > | may be considered, where 7 = 1 cor- 
responds to elastic case. This relation shows that np is the rigidity of the correspond- 
ing elastic material when 1 is the rigidity of sandy material and hence if » is considered 
to be the rigidity of the elastic material then the rigidity of the corresponding sandy 
material will be »/n. Also the generalised Lame’s constant may be defined as 


A= Ealy (1 +c) (1 — 20), 


With the above assumptions and taking into account of the principle of incre- 
mental deformation given by Biot. The Stresses and displacements produced in a 
dry sandy half-space under normal initial stresses due to a normal moving load on the 
rough surface of the half-space have been derived in this paper. The subsonic, tran- 
sonic and supersonic all the three cases have been discussed. The conditions for insta- 
bility, due to high stress concentration and the lines of generation of cracks have been 
obtained. The results in traasonic and Supersonic cases have been obtained in terms 
of Heisenberg delta function and are shown to coincide with the classical result (ny =1, 
S\1 = Sss = 0) obtained in Cole and Huth® in subsonic case and with Ghosh® (when 
h — eo) in transonic and supersonic cases. In subsonic case the numerical results for 
stresses developed have been calculated and represented by graphs for different values 
of sandy parameters 7 and the initial stress parameters J; and /,, 


GOVERNING EQUATIONS AND RELEVANT SOLUTIONS 


We consider a homogeneous, isotropic and linearly elastic sandy half space under 
initial stresses S,, and S;3 along x and z directions respectively. The half Space is sub- 
jected to a normal load F, independent of » and moving with a constant velocity v in 
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the positive x-direction. The moving load induces a state of plane strain in the half 
space whereby the y component of displacement vanishes and the remaining displace- 


ments and stresses are functions of x, z andtonly. The surface of the half space is 
assumed to be rough. 


i In the absence of any body force the dynamical equations of motion under in- 
itial stress P = S3; — S,, for two dimensional problems may be written as” 











0S), Css; OW y Cy 
ox = Oz a 4 ez. p ate (a) 
281 A833 Ow, ew 
= sata P ae ae (b) 


The stress-strain relations for sandy medium under initial stresses may be taken 


as 
agen A Qh eu A aw } 
Bene ren Bie AS ou | 
mo(2+st)P +o f 
| 
= of ied ou ow | 
Ig, = Sis = : ( az a ax ) ; seek A) 


where s;; are incremental stress components, ,, the rotational components, ¢1; the 
strain components and 7 is the sandy parameter A, }: are L_ame’s constants for the elastic 
material. 


Now from (a) and (b) and (1) the equations of motion for sandy medium under 
the considered initial stresses can be written in terms of displacement components as 


Cru 


ice Cu 
(A+ 2h +P) a3 ) 


+ (e+ na az? 


nP 
; +(A+o + 1 | 


aw Ou 
a 02 
dx oz Pon 2 


7 P \ aw a°w ( 1p 
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e2u o°w 
*2 ao Ot eel A 
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The boundary conditions at the free surface (z = 0) prescribe the normal stress to be 
delta function and the tangential stress to be balanced by the frictional force 


ie Afs = — F8(x — vt) 
Af, = — RF8(x — vt) A? 
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where, Afz, Afx being the incremental normal and shear boundary forces for unit init- 
ial area given by”? 





1) 6 A 71 

yeas (A + 2p) + ( +515 ) 3 

n oZ n Ox \ 

=oy (ey Bc (+ Pit eee) Come 
Afs = ( +5 ) dz a 2 5) ) ax J tS) 


Further, R being the coefficient of static friction at the surface and 5 (x) isthe Dirac 
delta function of argument x and is defined by 


ioe) 
On | ettx dk. 


—0oo 


s(x) = 


We take the steady state solutions of eqns. (2) and (3) in the form 


oo 
us j [A cos k (x — vt) + Bsin k (x — vt)] e-*# dk 


w = f [Cocos k (x — vt) + Dsin k (x — vt)] e-*? dk. Wet6) 
0 


-——+-—- 


Putting the expressions (6) for u and w in eqns (2) and (3) we get four equations con- 
necting A, D and B, C which are consistent only if 


gt +{ 2 > ou + oP) 4 ne vt — (u —$ nP) 
(4 + $ »P) (A + 2) 


(A + »w + 3 P)? 


TG FFP) O+2y) 


feales [Py 2= 21 A+2n+ yP)] npr? — y (vu — 4 yP)) 0 
G+EWP)AF RO 
If gi and q3 be the roots of equation (7), then uw and w can be written as 


qz 


co — z 
mei [{A,e ie + Ase ms } cos k (x — vt) +{B, : 


—kq z 


+ Be *} sink (x — vt)) dk SA ty 


fora) —kq z —k 
vid [{m, Bhe ‘+m, Bie "!'} cos k (x — vt) 


—kq z 


—kgz 
—{m Ave * + m,A,e “*} sin k(x — (vt)] dk (9) 
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where 


(A + 2u + aP) — Pov — (u + } nP) G3, 
(A +» + 4 nP) ai, 3 


m,: = 





Equation (7) is of the form x? + bx + ¢ = 0. It is observed by calculation that 
b* — 4c > O and in order that qi and q? may both be positive, the expression — b 
+b: — 4c > 0, which leads to the condition that 7p v? should be less than both [A + 

P 

2u + »P] (= T,, say) and [ — = |(= 7, say) and [{Pvr'n — (A+ Qu + Po)} 
x {v'4P — (4 — $ nP)}] (= c, say) should be positive. If g7 and q3 are both negative 
then Pf v’ » should be greater than both T, and 7, and also c is greater than zero. In 
this case q, and q, both are imaginary. In the isotropic case (1 = 1, S;, = Sss = 0), 
in order that g? and q? may be real, the velocity of moving load should be less than 
shear wave 6 or should be greater than «. When P v°x lies between 7, and 7,, g? q? 
becomes negative. So either g? or g? becomes negative in that case, i.e. either ¢, or 
g2 is imaginary. The above three cases are known as subsonic, supersonic and transonic 
case respectively. We shall discuss all the three cases in this paper. 


Inserting (8) and (9) in boundary conditions (4, 5) we get the four equations for 
the constants A, and 4,, B, and B, from which the constants are obtained as 


ye 2FRn [ (A + 2p) os) 
ee ices ee 


A 
re 2a [ (A+ 2u) mgr — (A + S33) ] 
iS imo Anke A* 
Bow — FA_[ (Qu — Sas — Su) my + (Qu + P) Qs ] 
SDS. Py a DP epee oa ake | cerca o 


...(10) 


oat (Cp. — Sys — S11) m, + (2p +P) | 
ak A 


where 
A* = [(2u — Sys — Sy) m + (Qu +P) qi] (A + 22) me gg — (A + Ss;)] 


— [(2u — Ss3 — S,,) m, + (2h + P) ge) ((a + 2p) m, qi 
— (A + Sss)]. (1) 


ANALYSIS 
: ] iP 
Case I Subsonic case (p vi< elke G4 
In this case P v*y is less than both 7, and 7, and both the roots gq, and q, are 
real, Inthe isotropic case (y = 1, S;,; = S33 = 0) this corresponds to the case that 
the source moves with a velocily less than velocity of shear waves. 
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Thus, from (1) using (8), (9), (10) and (11), the expressions for the stresses Sis 
and S33 may be written as 











apa mo { Gs io ae box — vt) 
ae | ae f3bm \| (12) 
and 
Sas = re eee PeG te? + 2R | oe 
- rast ts » | (13) 
where 


M, =m + qi, M, = m, + q, 
Or = (x — vt) + gf 2,0,=( - +2 
G, = P3m + Py qi, G, = Py m, + Py Gz 


Py = (A + 2p) m qi — (A + Sasa), P2 = (A + (24) me gy — (A + Sa3) 


i 


P; = (2 — S33 — Sn), Ps = (24 + P), Ps = (A + 2p) mq — A, 
P, = (A +- 2p) mez qe -—* ‘AS 


A* = 0 is the frequency equation of the Rayleigh waves in the medium under consider- 
ed initial stresses. Hence, when the velocity of the moving load coincides with the 
velocity of Rayleigh waves obtained from A* = 0, the stresses will be infinitely large 
and fracture will take place in the medium. 


In the absence of initial stresses (i. e. S33; =S,, = 0) when 7 = 1, the developed 
incremental stresses given by (12) and (13) reduce for smooth surface (i. e. R — 0) and 
the results are seen to coincide with the results obtained by Cole and Huth*. The 
Stresses at X = vfi.e., at the point directly below the load may be obtained as 








see — FR [tm + @) (0+ ma — (06+ 1} 
n A qe 
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4" 
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where 
= Af2u, 1, = $1:/2p, 1,“= S33/2p. 


Case I1—Transonic case :[i.e. when <(u -- “ ) << p v2 


] (= T,) 
Ae la a tae arome ] 


Next, we consider the case when Pv’, lies between T, and T,. In the isotropic 


case this corresponds to the case that the source moves with a velocity greater than the 
velocity of shear waves but less than the velocity of longitudinal waves. 
In this case 7, is imaginary and q; is real. Here, we replace q, by iq, where 4) 


is real. The developed incremental stresses S,3; and S;, are given by 








Fp x —vt M, G, 
aL G, “ers Y {2m 3, {(x — vt) — gq; 2} 


— 2n 5, {— (x — vt) — q; Z}} 


5)3 = 


9 o 
+ 2R a CAE heat [27 5, {(x—vt) — J, Z} 





+27 {—@—w)— 4, 2p | (16) 
ii | Gate Fate Pe St tm 8, ((x — vt) — gf 2} 
+ 2m be(-@—w)- giz} b+ arf BPE 
% are [2a 54 {(x — vt) —qi 2} 
pata) 


— 2m 8, {(x — vt) — qj 23] I 


where, 5, (x) is the Heisenberg delta function and is given by 
5, (x) = d(x), x =0 


—_= — ee ee 0: 
27 ix 


The Heisenberg delta functions within the brackets may be written in explicit form as 


[2a 5, {(x — vt) — gy 2} — 2m oy {— (x — vt) — q) Z}] 


(equation continued on p. 172) 
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and from equation (7) g; may be written as 








, -{ b+V/b? — 4c fe aitts5 
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1 
From expression (17) with the help of (18), it is observed that when —~( » — 


ut ) <py? <— (A + 2» + Pn), S33 becomes infinite along the lines (x —yf) 


=-+4q, z. Hence crakcs are produced along these lines. From (19) it is seen that the 


slope of the crack lines q, depends on the initial stress parameters /, and /,. However, 


the initial stresses have no effect on the cracklines when J, = I, i. e. the initial stresses 
are hydrostatic in nature. It is easy to calculate the slope of the cracklines from equa- 
tion (19) for particular values of the initial stress parameters and at any time ¢. 


Case I1]—Supersonic case : py > (A + 2 + Pn) Tt) 
Consider the case when p v*y is greater than T,. 


= I, = 0) this corresponds to the velocity of source is 
well as P-wave velocity in the medium. 


In the classical case (w= 1,7, 
greater than the shear wave as 


In this case g; and 92 both are imaginary, Replacing q, and q2 by iq; and ig; 
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respectiv * , i 
pectively, where q; and q; are real, the expressions for 5,3 and s,, are obtained as 


Fu M, G, 
At i {27 5, (x — vt) — gq, z} — 2m 8, {— (x — vt) 


Ssi3 = 











— qi 2h} — “ES! (an 8, (x — vt) — gi 2} — 2 8g (— (x—Wf) 
— 9, y+ 2R {PM Dow 8, (em) - 442 

+ 2n 3, (x — vt) — 9) 2) — “2 [oe 8, (x — 1) — a 2} 
eet ix va it ] ...(20) 


and 


F P,G 
S33 = Ke ae {20 8, {((x — vt) — gy z} + 278, {— (x — vt) 


—q;, z}}— Fs St (29 8, {((x — vt) — q) z} + 2m 8, {—(x—vt) 


; Pak: 
— gf.2)) + 2R {FP fon 8, (x — vt) — af 2-208 
re 
(—@ —m)— 9, — FEPeL om — 9) - 
— 2n 8, {— (x — vt) — 9; ait] 148) 
where g, as obtained from (7) is given by 
; b —vVb? 4c 

From eqn. (21) it is observed that when the point source moves with a velocity 


+ a+ 41?) 
greater than | 4 p 





pat then four cracks are produced along the lines 


(x — vt) =~ +9102 instead of two cracks as obtained in case II. Beyond the crack- 


lines the stresses may be calculated from equation (21) using the definition of Heisen- 


berg delta function. These cracklines will be very much affected by the presence of 


initial stresses because their slopes + q},, are functions of initial stresses present in 
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medium. The slopes of cracklines may be calculated from (18) and (22) for different 
values of initial stress parameters and elastic constants. 


Values of 


NUMERICAL CALCULATION AND DISCUSSIONS 


The values of s33/F and s,,/FR have been calculated (ICL 1901A) for subsonic 
case when z = 100 from eqns. (14) and (15) (in dimensionless form) in elastic (nj = 1.0) 
and sandy materials (jy = 1.5) taking A = y, as a particular case, due to increasing 
velocity of moving load in presence of pre-comperssive stress along x-direction (/, =0.0, 
— 0.2, -- 0.4, — 0.6) and pre-tensile stress along z-direction (/, = 0.0, 0.2, 0,4, 0.6, 
0.8) and also for the case free from pre-stresses to facilitate the comparison. These 
results have been presented in Tables I and II. It is observed that the stresses produced 
in elastic medium is always less than that in sandy soil. 


TABLE I 


normal stresses and shear stresses due to increasing velocity of moving load in an elastic 
as well as sandy half space under compressive stress along x-direction. 


z=100, ,= 
I; v/B $ao/F $1:/FR 
= 1.0 n= 1.5 = 1.0 7 = 1.5 

(1) (2) (3) (4) (5) (6) 

0.0 0.1 0.641 x10-2 0.643 x 107? 0.242 x 1074 0.366 1074 

0.0 0.2 0.655X10°2 0.555 x 1072 0 1011078 0.15610-3 

0.0 0.3 0.681 X 1072 0.707107 0.246 1073 0.396 10-3 

0.0 0.4 0.723 x 107? 0.781 10-2 0.490 10-8 0.850 1073 

0.0 0.5 0.790 10-2 0.920107? 0.903 x 10-3 0.178 x 10-2 

0.0 0.6 0.902 x 10-2 0.124X1071 0.16510°2 0.421 1072 

0.0 0.7 0.1111071 0.275 x 1071 0.321 x1072 0.177>'1071 

0.0 0.8 0.167107! infinitely large 0.778 x 1072 infinitely large 

0.0 0.9 0.761107! infinitely large 0.650 1071 infinitely large 
—0.2 0.1 0.833 x 107? 0.101107! 0.122«1072 0.220 1072 
—0,2 0.2 0.862 1072 0.108107 0.1391072 0.265 x 10-2 
—0.2 0.3 0.918 x1072 0.124x1071 0.172107? 0.365107" 
—0.2 0.4 0.101107} 0.162107! 0.232 1072 0.61610~-2 
—0.2 0.5 0.119 1072 0.307 x 1071 0.346 1072 0.165107? 
—0.2 0.6 0.156107! infinite'y large 0.608 x 10-2 infinitely large 
—-0.2 0.7 0.284107! infinitely large 0.158107 infinitely large 
—0.4 0.1 0.127X1071 0.361 x 107! 0.368 1072 0.171107 
—0.4 0.2 0.1361071 0.578 1072 0.420 10-2 0.305107! 
—0.4 0.3 0.154107 infinitely large 0.534 1072 infinitely large 
—0.4 0.4 0.194107 —do- 0.789 x 10-2 —do- 
—0.4 0.5 0.3111071 —do- 0.158107 -~do-— 
—0.4 0.6 0.209 —do-— 0.207 =—do— 
—0.6 0.1 0.347107 —do— 0.1632071 do 
—0.6 0.2 0.454X1071 —do— 0.228 x 107) —do— 
—0.6 0.3 0.192 


hoe 0.579 10-1 a Pe 
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TABLE II 
Values of normal stresses and shear stresses due to increasing velocity of moving load in an elastic 
as well as sandy half space under tensile stress along z-direction 
z= 100, J; = 0.0 
SS 








I v/B $32/F $i3/FR 

n= 1.0 7 =1.5 7 = 1.0 n= 1.5 
(1) (2) (3) (4) (5) (6) 
0.2 0.1 0.653 X10 2 0.686 107" 0.244 1074 0.3911074 
0.2 0.2 0.670107 0.714107 0.1031078 0.170 10-8 
0.2 0.3 0.698 X10 -- 0.76910 0.251 1073 0.445107 
0.2 0.4 0.746 107° 0.874107> 0.510 1073 0.102107? 
0.2 0.5 0.824107? 0.110107! 0.968 1073 0.244 X10~? 
0.2 0.6 0.962 107? 0.200107! 0.186107? 0.900107? 
0.2 0.7 0.126X1071 infinitely large 0.407102 infinitely large 
0.2 0.8 0.242K1071 —do-- 0.148101 —do— 
0.4 0.1 0.732 1072 0.9181072 0.284107! 0.601 x 10-4 
0.4 0.2 0.754107 0.98310 2 0.1211073 0.276 1073 
0.4 0.3 0.795 X10-2 0.113K10~! 0.30310-3 0.816 1073 
0.4 0.4 0.866107? 0.1501071 0.640 10-2 025310> 
0.4 0.5 0.994107 0.512107! 0.131107 0.351107} 
0.4 0.6 0.126107! infinitely large 0.295102 infinitely large 
0.4 0.7 0.222*1071 —do— 0.107107! infinitely large 
0.6 0.1 0.912X10°7 0.219X10~4 0.3951074 0.269 10-3 
0.6 0.2 0.952107 0.323 X1072 0.173 1073 0.232X10-* 
0.6 0.3 0.103107} infinitely large 0.458 x 1073 infinitely large 
0.6 0.4 0.119X1071 —do-- 0.108 10-7 —do— 
0.6 0.5 0.155107} —do-- 0.28310" —do— 
0.6 0.6 0.359X107! —do— 0.186107} —do— 
0.8 0.1 0.139107 —do— 0.789 10-4 —do— 
0.8 0.2 0.1§2X10-2 —do— 0.378 X 1073 —do— 
0.8 0.3 0.185107! —do- 0.126X10~ —do-—- 
0.8 0.4 0317%107 —do— 0.643X10°2 —do— 


a 


TABLE III 


Critical values of v/® at which A* = 0 
nn LEE 





I I; v/B oak 
ea 0 iVrsebe! WSS 
—0.4 0.0 0.6098 0.2920 
—0.2 0.0 0.7838 0.5721 
0.0 0.0 0.3187 0.7505 
0.0 0.2 0.8625 0.6650 


7 See eee 
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Fic. 1. Development of normal stress 533 in elastic half space (y = 1.0) under compressive pre- 
stresses along the x-direction. 


The critical values of v/8 have also been calculated from A* = 0 at which the 
stresses developed will be infinitely large creating fracture in the medium. The Critical 
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Fic. 2. A comparison in the development of normal stresses in pre-stressed and free from pre- 
stresses half space due to increasing velocity of moving load. 


velocities of moving load giving rise to fracture fora few values of pre-stress para- 
meters in elastic and sandy material have been presented in Table III. It is found that 
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Fic. 3. Development or normal stresses in an elastic half space moderate and high tensile-stress 


z-direction, 
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Fic. 4. Critical velocities of the moving load creating fracture in a pre-stressed elastic half space. 


the sandy material gets fractured at lesser velocity of moving load in comparison with 
the elastic material i. e. sandy materials are less stable to moving load than the elaStic 


one. 


Figures 1-3 to show the nature of the curves of the normal stresses (s33;) develop- 
ed in a pre-stressed elastic halfspace due to increasing velocities of moving load for 
compressive stress along x-directions and tensile stress along z-direction. It may be 
that the curves are asymptotic for certain values of v/3 depending upon the pre-stress 
parameters /, and /, suggesting that the stresses are very high creating instability in 
the medium causing fracture at the corresponding critical velocity of the moving load. 


Observing Figs. 1 and one infers that high pre-comp ressive stresses along 
x-direction will cause the instability at lesser velocity of moving load than the classical 
case. From Figs. 2 and 3 the same remarks can be made for the case when the elastic 
medium is under pre-tensile stress along z-direction. Also from the above figures it is 
clear that as the magnitude of compressive stress along x-direction or tensile stress 
along z-direction increases the fracture takes place at less and velocities of the moving 


load. 


Figure 4 represents the curves for the critical values of v/8 at which the stresses 


developed ina pre-stressed elastic material will be infinitely large creating fracture in 


the medium. 
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THERMO-ELASTIC WAVES FROM SUDDENLY PUNCHED 
HOLE IN STRETCHED ELASTIC PLATE 
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Coupled thermo-elastic wave problem of infinitely extended elastic plate of 
thickness d (subject to initially a state of axially symmetric hydrostatic ten- 
sion) a flate nose cylindrical Projectile of radius C travelling with velocity V 
strikes the plate and begins to punch out a hole of radius equal to its own, has 
been studied under cretain assumption. The integral transform technique is 
used. The expression for stresses and temperature are derived. Long time 
and short time solution are presented. 


INTRODUCTION 


The elastic deformation due to time dependent surface traction (from Duhamel— 
Neuman analysis) can be determined without reference to its thermal state. This pro- 
cedure is adopted in the solution of dynamical problem in the classical theory of 
elasticity. The elastic constants which appear in the equation of motion are defined 
under adiabatic condition and since elastic waves being non-dispersive produce no 
increase in thé entropy of the solid, no inconsistency seem to arise. This description of 
elastic wave propagation is physically over simplified. In fact a change in volume must 
cause the temperature as well as stresses. When a longitudinal wave passes through a 
solid the elements are successively compressed and dialated. These phenomena are 
accompanied respectively by heating and cooling. Since the thermal conductivity of 
the solid is non-zero and the disturbance has finite frequency, the source of energy 
will be converted into heat energy during the first half of an oscillation will not be 
recovered during the dialation phase. 

With the foregoing remarks in mind this paper is prepared. A formal solution of 
coupled thermo-elastic problem of infinitely extended elastic plate of thickness d sub- 
ject to initially a state of axially symmetric hydrostatic tension i.e., 6, = og = A,a flate 
nose cylindrical projectile of radius C travelling with velocity V strikes the plate and 
begins to punch out a hole of radius equal to its own has been obtained. The problem 
is solved under certain assumptions : 

(a) The plastic flow due to punching 1s very localized to the neighbourhood of 


the punch sections. 
(experiment also support for V > 2000 fps). 
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(b) The punching begins instaneously at ¢ = 0 over the whole punched section, 
based on a small value of plate thickness d and large value of impact speed V. 


(c) The punching action takes place at velocity V/2, the particle’s velocity in the 
compressional wave that develop in both projectile and plate on impact i.e., the plate 
material below the projectile is removed as a plug at V/2. The corresponding punching 
time is therefore 2d/V = /, based on large ratio of diameter of projectile to plate 
thickness. 


A few problems are solved on thermo-mechanical coupling effect! Kumar un- 
coupled problem Miklowitz*. The result of Miklowitz® are obtained as particular case. 


2. MATHEMATICAL FORMULATION AND SOLUTION OF THE PROBLEM 


The origin of the cylindrical co-ordinates (r, @, Z) is taken on the axis of the 
cylindrical hole. For radially symmetric case the only non-zero displacement compo- 
nent is wu (r, f) in the radial direction for plane stress problem. 


The thermo-elastic equation of motion in the absence of body forces can be 
written’, 
O Grr ] (Sr, = Soe) Z: au 


ot r — Poe (1) 





and conduction equation 


oT de 
2 7T ae 
PC,K y* T= PC, = +m T, ‘ em 


where 


eet Nae oF 
Vv peta eats > ap. = (3A + 2p) %. 


The strain-displacement relation and stresses are : 


u u 
fe ee te ee 
lg 


or 


; : 
Or = Ae + 2p = — (3A + 2p) « (T —- 7,) 


Seo = Ae + 2p — — (3A + 2y) a, (T — T) 


5 Cu u 
€ = rr + Cog = or. + oF 5 GD 
where A, » are Lame’s constants; «, is the coefficient of linear expansion; (T — T,) the 
deviation from the equilibrium temperature 7(; P the density; C, the specific heat at 
constant cubical dialation; and K the conductivity of the material. 
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The boundary conditions of the problem are: 





























o;,,r = 0; < 0 : 
= - Ho<r<i { oe 
| 2d 
er te l= —. ...(4) 
Attic j V 
From the relations (1), (2), (3), we can write equation of motion as, 
A+ 2n @u A+ 2h @ u (3A + 2p) 7 
Par? Picov ary pi p a a (T — To) 
au 
Sia gas) 
and Conduction equation : 
2p a PCy OF mIo = ou =) 
ie aiey aie K ot\ or r ) at 
By Introducing the following non-dimensional quantities 
GC PC? u 
Ci > 1 > 1 
ae te ie Ty a 
T, = T, Ct = ar ee as. G7) 
c= a is the coefficient of thermal diffusivity (7). Equation (5) and (6) reduces to 
avU 1 aU U o Ure 2.87; 
aR? | ROR RF” Tee? "OR ma 
a T, [eeT i ery soa =) 9 
SARE ORG SOx7 BAKER 2 R re 
where 
_ 403A + 2m) as}? Ty is a coupling factor. 
eet (1S Dad hehe oe te 
Equation (9) and (8) can be written as 
a (aU U ) eU ) 
of Berek — —T,)= : 
(a) on (OR cake to aE ‘ ...(10) 
F I oT, wv) a Mas a 
(f+ BB) 


We define transformation 
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7} 
= oo 
(a) f (R, p) = J f (R, +) e-?* de | 
| 
and its inversion is : hank 
; c+1o0 
(b) f (R, 1) = = | F(R P) ev" dp. | 
c—100 ' 
J 
Taking transform of eqns. (10) 
Bi (BUM, VO ners i } 
(a) axl ae + R -7,) = p> U | 
» (4 ae Ce o rie J f --(12) 
() (ore) or P a ees j 


where ‘Bar’ stands for transformed function and pis transform variable with the 
assumption of 
é a 1 
ge ee, «ee sat . as 1 
D aR: D ( aR + z) enw as DD 
D, = Di D, 


Equations (12) can be written as 


(a) (DD' — p*) U = DT, 
(b) (D' D — p)T, = 38 pD' U. pe 
Operating D,, D, on (a) and (b) respectively, we get 
| Bt - (mi + mt) Dy + mi m: | O=0 ...(14) 
| Of —( mt + mt) De + mi m; | 7, =0 s3C15) 
where m; and mi are the roots of the equation 
m* — m*[p* + (1 + 8) p] + p? = 0, (16) 
Equation (14) is equivalent to 
(a) (D, — mt ) 0, = 0 } 
| 
\ 
i ee . (12) 
(b) (p= m ) 0, = 0; Wate | 
J 


and eqn. (15) is equivalent to 
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(a) (D, — m: ) 7, = 0 } 
| 

\ ..-(18) 
[ 
J 


The general solutions consistent with the boundary conditions are 


(b) (0. = m; iss =O; T, = 7, + 75. 


U = U, 4. U, = AK, (m R) + BK, (m, R) } (19) 


T, = Ty, + Ti, = EK, (m, R) + FKy (rm: R) 


where A, B, E, F are arbitrary constants and Ky (Z), K, (Z) are modified Bessel func- 
tion of second kind of order zero and one respectively. 


Substituting the values of UT, from (19) in eqns. 13 we get the following rela- 
tions in constants 


es See ees ee ey ...(20) 


and thus 


U = AK, (m, R) + BK; (m, R) 


(p= mi ) p? — m: 

T, = ee A K,(m, R) + a B K, (m, R). 
As the displacements are known, we can compute the components of stress and 
temperature. 

orr = A Mi (m, R) + BM, (m, R) 

T, =AN,(m, R) + BN, (m, R) ae § 
where 

B? = (u/A) + 2p 


2 
M, (m; R) = 2 Ko (mm R) + BK, (m, R) 


M, (m,z R) = - Ky (m, R) + B? Ki (m, R) 


(p? — m; ) 
Rui RY = ——— Poh R) 
(p? — m, ) 


aos 2ap 
N,(m, R) = ear Te Kg (mz R). ( 


2 
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The rest constanss A, B can be determined by using the boundary conditions (4) after 
taking their transformation, and equating (21), we get 


Bol 
A M, (m, b) + BM, (m, b) = Pre (e'? — 1) 


A N; (m, 6) + B Nz (m, b) = 0. Ae EY 
Solving above two equations for A and B, we get 


A (lL — e-!?) Ne. (m, b) 


A = “Ip [M, (m, b) N, (m b) — M, (m b) N, (m, by] 
— — A(1l — e7'?) N, (m, b) -..(24) 


Ip* [M, (mz, b) Ni (mm, 6) — M, (m b) N, (m, 6)) 


whereb = “1 a 
x 
After substituting the values of A and B in (19) we get displacement and tem- 
perature in image space 


A(l — e-!?) 


O = — 
I p® F (p) 


[M2 (m2 b) Ky (m, R) ~— N, (m, b K, (m, R)) 


B ey eae) 
Ty Ae oe ae 





I p® F(p) m, 
Po 52) 
aoe m N, (m, b) Ky (m, R) | 
2 
where 
F(p) = M, N, — M, Nz. t2o) 


Long Time Solution 


. 


The long time solution can be obtained by expanding the roots m; ,m, of (16) 
for small value of p. The expansion is done in Taylor’s series* 
m = (1 + 84) Vp + O (p)3 
m, = (1 + 8) -'/? p+ O (p)? BALES), 


Neglecting higher power of p and denoting as 


m = 1, Vp where = (1 + 8)712 
Lp k= (1 + 8)7 72, (27) 


Me 
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By replacing the values of m, and m, from (27) in (25), we get displacement and tem- 
perature in image space for long time solutions, 


Coe ae TN hp BY Kt Ve RYN hd) Ree 
I p® F(p) 2 (‘oP iG, VVv> — N, (1, p 6) K, (; p R)) 
- — Al — ee) (p? — Ii p) : 
Ty = ~ 1p? F(p) Ep NG Po) Koh V pR)- 
(pie feet} 2 
on N, (lL, Vp b) Ky (I, p R) | hee) 


and similarly the stress components or (R, P), coo (R, p) and temperature 7; (R, p) 
can be evaluated in image space. 


Short Time solution 


2 
The short time solution can also be obtained by expanding the roots Ua) @ 


(16) for large value of p. The expansion is done by Laurent series‘ : 


s(1 — 38 s(1 — 36 + & 
mi = p++ + C=) 4 “UU BEN +00 
. 5(1 — 3) AN Eee ee ae IS, Oy, s 
en ae rapa ee op nO 
or 
8 8 . 
m—=p te + C=) 4 or | 
> ...(29) 
= eal - 3/2 ve 
NP ON parr) j 
Neglecting higher powers of p,m, = /,; m, = /, where 
8 $(4 — 8) 
GAIT toe ae le Nema Ta 
ue ...(30) 


l= VP ae Qa/p° 


By substituting the values of m, and m, from (30) we get corresponding displacements 
and temperature for short time solutions in image space. 





7 A(l —e) ry 7 by KL (hh R) — Ni (hh) Ki (a RD] ...(31) 
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A (1 — ee”) F (p? — I’) R 
1= | p* F(p) aoe OR ah ) 


I 
| 


(p?— I) 
= pr Ni a R) | ...(32) 


and similarly stress components, orp, oe¢ can be evaluated in image space. 


The stress components ore (R, t) and T, (R, t) can be obtained by taking the 
inversion of the transform. 


' c+too 
Orr (Ri t) = a | Grr (R p)eP* dp. 
¢c—ico 


Similarly all other components and temperature can be evaluated. This formally com- 
pletes the problem. 


3. DisCUSSION OF RESULTS 


We have discussed the above problem under thermo-mechanical coupling effects. 


In eqns. (17) we have noticed that there appears one more displacement com- 
pared to uncoupled problem? which is responsible for thermal waves in the medium. 
As a matter of fact consideration of Thermo-mechanical coupling is desirable for all 


elasto-dynamic problem at least for high speed impact problems otherwise we neglect 
a good amount of energy. 


The thermo-mechanical coupling is done by a coupling factor 5 and if 8 = 0 the 
temperature from (19) will vanish. 


If 5 = 0 the displacement function agrees with Miklowitz’? in image space and 
thus con-sequently stresses. 


If’ = 0,t = Othe problem reduces to uncoupled stationary problem. 


- 
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NOTE ON MINMAX PRINCIPLE FOR HEAT CONVECTION EQUATION 


M. A. GoPALAN 
Mathematics Department, National College, Tiruchirapalli 620001 


(Received 24 February 1988) 


A variational principle for the energy equation with viscous dissipation sub- 
ject to Dirichlet condition is obtained. 


1. INTRODUCTION 


Barret et al.' have given a minmax principle for a non-self adjoint system and 
applied it to Navier-Stokes equations. 


Here, we give a minmax principle for the heat convection equation with viscous 
dissipation when the surface temperature is prescribed. For simplicity, we mention 
the situations where we may obtain the saddle functional. 


2. MINMAX PRINCIPLE FOR ENERGY EQUATION 


The differential equation for the steady distribution of temperature in a viscous 
incompressible liquid is given by 


m I 
VT =Pc,V.yT — —o EC 


where « is the Thermal conductivity of fluid, 
C,, the heat capacity per unit volume of fluid, 
J the mechanical equivalent of heat, 
® the viscous dissipation function, and 
V the fully developed Jaminar velocity of fluid. 
The boundary condition considered is 
T = T, on the boundary S. Vee eed 


Introduce the functional 
J (Ts, Ts) = Sie (v T,)* — (9 7)) — 72.9 7; 
y 


+ 7,V.9 Tz, — 2a (T, — T:) 0] dV voe(2.3) 


190 M. A. GOPALAN 


Subject to 
T, = T, = T,onS .--(2.4) 
where 
« = («JP C,), « = (1/P C, J). 
Equating the first variation in J to zero, we get 
¥7T,=V.yT%—«® PAE ESY) 
<V7,=V.9gT, —<«9. ...(2.6) 
It is easily seen that 


T, =T,. 


Then (2.5) and (2.6) reduce to (2.1) and finding the stationary point of the functoinal 
(2.3) is equivalent to solving the original problem. 


If 7, is fixed, the second variation in J is 


82 J = f «(yet aV (27) 


where & is an admissible function. Keeping 7, fixed, we have 


SJ = — J (vy n)? dV (2.8) 


where » is an admissible function. From (2.7) and (2.8) we notice that if T, is fixed, J 
is minimum with respect to 7, and if 7, is fixed, J is maximum with respect to 7. 
Therefore, J represents a saddle functional and the solution is obtained from 


Max Min 
i, I, /(7,,T,) ...(2.9) 


where the trial functions satisfy the appropirate boundary conditions. The following 
solution algorithm is suggested : 


(1) Choose T, to satisfy T, = T, on S. 
(2) Solve 
Ovi Tl =v his a ®, 7, = T,onS. 
At the end of this step, the corresponding value of J is 


J (T; (T2), T,) = fe(vni-y T,)? dV ...(2.10) 
providing a direct estimate of the error. 


(3) Maximize J with respeet to T, and thus obtain the unknown parameters in 


De 
< 
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Note that for the time-dependent temperature equation 
6T 
WT = =, + V.yT —a® atin) 


with the initial and boundary conditions 





T=Oatrt=0 
the saddle functional may be found to be 
J (T,, T;) =JIk(y L(y Tey) —27, V.VT; 
+7,.0. 9T, — 20 (T, —T)0— T, an 
oT ae | av at. tek 


3. APPLICATIONS 


(A) Generalized Couette Flow with Suction and Injection 


Consider the steady two-dimensional flow of a viscous incompressible fluid bet- 
ween two parallel flat plates. The equation for the temperature distribution with 
viscous dissipation is* 


d? T* ar® aut 
Fm = PR Tre PE ( S- ) ih Sel) 


subject to the conditions 
PO) = 0/T* (1) = 1. yal oak) 
The saddle functional is found to be 


1 


aT, \? dT, \* aT, 
J(T;, T.) =|(( ae ) ag ( dn ) — PRT; dy 
0 





d: ae 2 
+ PRT, — 2PE(T: — 7.) ( a ) dn, (3.3) 


The algorithm similar to section 2 is used to obtain the solution and the estimate of 


the error is given by 


1 





dT, dT, \? 
5(7, (7) T) = — \( F2 - Ge) an 3.4) 
0 
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(B) Parabolic Flow Between Two Semi-infinite Parallel Plates 


Consider the parabolic flow between two semi-infinite parallel plates y = 0, 
y = 2h maintained at a constant temperature To, the temperature of the incident fluid 
being zero. Neglecting the axial heat conduction and viscous dissipation, the governing 
equation for the fluid temperature is® 


Oa oT 
K Oy = U (y) Ox cuataeo) 


where U (y) represents the fully developed laminar velocity in channel. The saddle 
functional is seen to be 


h 


serwty = [{[(#Y-( By) - on B 


+ U(y)T; ace ] dy dx ...(3.6) 





and the corresponding error estimate is 


00 gh 
I(T: (T), (T) = — \| ; (So isc Os 
PF dy oy 





)a dx. ata} 


We notice that for the slug flow between parallel plates, the saddle functional is 
represented by (3.6), where 


U (y) = VU, (constant). 
(C) Magnetohydrodynamic Pipe Flow 


Consider the steady flow of a conducting liquid along a pipe under the influence 
of a transverse magnetic field. If S represents a cross-section of the pipe in the (x, y) 
plane, the fluid velocity w and induced magnetic field h satisfy‘ 


dh 


2 ——_—=—— i ic—-—— i 
Vere ies lin , ...(3.8) 
Ow 
*h+M—=0i 
Vih+ By OinS +06(3.9) 


with 
w=h=QondSs 
for stationary insulated walls, being the Hartmann number. 


Equations (3.8) and (3.9) may be uncoupled by addition and subtraction to give 


‘ oP 
UEP Aes cea 
oy l In s daslavll) 
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éQ : 
NO oR ar a a -lins (3.10) 
where 
P=w+A, Q=w-—h 
and 
P=Q=0on8S. Pom Bp PA 
The saddle functional for (3.10) is found to be 
J (P,P) = He 2 i 
(P,, ) = (v P,) (Vv P.)* + M P, oy 
OP. 
— MP, ay —2(P; = ») | ax dy (3@G43) 


and that the corresponding saddle functional for (3.11) namely, J (Q, Q,), is obtained 
by replacing M by (— M) P, by Q, and P, by Q, in (3.13). These principles are then 
used to obtain the approximate fluid velocity w and induced magnetic field / from 


w= 4(P + QO)andh = 3(P — Q). 
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Unsteady flow water behind weak and strong shock waves produced by the 
detonation of a spherical charge, is studied by the method of perturbations. 
Results are compared with those obtained by similarity methods elsewhere. 
Effect of gas bubble on the shock front is neglected in the present paper. 


1. INTRODUCTION 


The knowledge of flow profile behind the primary stock waves in water is of 
immense importance and also of academic interest. Using the method of similarity, 
Kochina and Melnikova'** have studied flow profile behind the shock waves, produced 
by the explosion and by the piston motion respectively, and by perturbation method, 
this problem is solved for flow behind shocks in air by Singh®. Assuming energy bet- 
ween the shock front and the piston surface to be variable, flow profile behind the 
shock in water has also been studied‘. 


In the present paper, we have studied the unsteady motion of water behind the 
shock waves, using perturbation method. Law of attenuation of shock was studied 
earlier in a series of papers by the first author®-8, theoretically as well as experimen- 
tally. In the above papers and also in the present paper, effects of gas bubble on the 
flow are ignored. 


Flow behind shock waves is governed by the equations of motion of compressible 
fluids which are integrated by the methods of perturbations, taking shock front as one 
of the boundary. Equations of motion are reduced to two differential equations in 
the non-dimensional fluid parameters f(A, A p), g (A, 4 P) and A, where f, g, and A are 
non-dimensional fluid velocity, density and distance respectively. Expressing para- 
meters f and g in the form of converging series, variation of these parameters with 
respect to A and shock strength Ap is obtained. Variation of parameters fand g 
behind spherical shock waves is shown in Figures | to 4 for the two cases of weak and 


strong shocks respectively. These results are similar to those obtained earlier* in 
piston problem. 


2. Basic FORMULATION OF THE PROBLEM 


We have earlier studied® the Propagation and attenuation of spherical shock 
waves, produced by the detonation of an explosive charge in water. It is our aim in the 
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present paper, to study the unsteady fluid motion behind the spherical shock waves. 
Basic equations governing this motion are, 


op a] 2Pu 





Or oT ee (P u) + a ec 0 eat) 

ou Gu 1 @ép 

Be ane’ pe ope (2) 
dE p dp : 

tdi abe dhs ...Q) 


Hugoniot equation of state of water is 
U=a+ bu v4) 
where a = 1.5565, b = 1.9107 (Madan®) and other symbols have usual meaning. 


If at any time ¢t, R is the radius of the shock front, the jump conditions across 
the shock front are’, 


P2 = Pia? d (8 — I)/{b — 8(b — 1)? at E), 
U = adl{b — 3 (b — 1)} ...(6) 
EY = ES +([(8 — 1) af{b - 8 (6 — 1)}p Ean tf 


where subscript 2 denotes values of fluid parameters behind the shock front and 
5 =P,/P, is the shock compression. Variation of shock compression 5 with the shock 
radius is given by’®. 


8 [a (8 — 1)/{6 — 8 (b — Ij = 3 « O If(4z P, R°) (6) 


where QO is the heat of explosion per unit volume, R = R/R,, Ry being the radius of 
undetonated explosive charge. Flow behind the shock wave is governed by equations 
of motion (!)—(4). Thus our aim in the present paper is to find the Solution of these 
equations with the help of boundary conditions (5)—(7). 


We define the shock strength by a parameter AP so that 
Ap (6, PP 8 — | 9) 


where AP is assumed to be very small for weak shocks and for strong shocks 0 < Ap 
< 1. Largest value of Ap is 0.7 for the available conventional explosives. Using ex- 
pression (9) in (8) we get after differentiation and expansion (Appendix A). 


CAP _ pag ...(10) 
> Talon BAP/R 


B = By +B, Ap + B, Ap? + ... 
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where 


B, = 3/2, 8:1 = — 3(2b — 1)/4, B, = 3 (4b — 1)/8. CU) 


3. DiscUSSION OF THE PROBLEM 
To evaluate the variation of fluid parameters behind the shock front, we de- 
fine non-dimensional parameters f, g and A as, 


ulU = f (A, AP), p/P, = g (A, AP), r/R = A sttes 


where f and g are functions of A and 4P, r the radial distance measured from the point 
of explosion and u,P the fluid velocity, density at distance are respectively. Substitut- 
ing the parameters from (12) in equations (1) and (2), one gets after some simplifica- 
tion (Appendix A) 


Saf "dj 8x = BAP sg SBAP TEAR) aie ke} 
ae ost BO — Yea 

(a oAp + fae) pap -(f- O- {Seay} 
(14) 





where fj, a, : &,, are partial derivatives of f and g with respectto A and A re- 


spectively. Since fand g are functions of \ and AP, we can write f and g in the form 
of converging series. 


f(A, AP) = fo (A) + fy (A) AP + fe (A) AP? + ... 3,015) 
& (A, AP) = Bo (A) + 8, (A) AP + g (A) Ap? + ... ...(16) 
where fo, 20, /;, 8; etc are functions of A only. 


After substituting the expressions (15) and (16) for f and g in (13) and (14) and 
comparing the coefficients of same powers of AP, one gets equations of different order 
as zeroth order equations. 


Be-Jeo te (So =A) Bir 2fngo A OU neha 

(fod) Sp tee aa O ...(18) 
First order equations 

8 fi + (fo — 4) 8; —A&, = 0 (19) 

8o(fo — A) f) +8; +4, = 0. ...(20) 
Second order equations 


Bot, +(fo—AS, —m = 0 mre 
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B(fo— A fi +e, +h, = 0 .t22) 
where 


h, =f, 8 fy + (26—2— 8, 871) 2" — By go(f, + bh) 
he = — 81 fy — fi 8 + Bo (Bo + &:) — 2( fi Bo tSoardla 
hs = f2 Bo fy theo, 
+ 8 (8) 8 — $285 — 8&1 (2b — 2) g,° + 26% — 6b 


+ 3] —g/ (g: 9; — 26 +2) .. (23) 


— Bo gol fo (b? — 26) + fi b +2 fy] 
— Bi Bo (b fo + fi) 


m=—£:f, —-8:f, —f2 8 —Sfi 8 


— Bo (89 — 8; — 282) + Bi (So + 21) 
= 2 (fo &2 + fie, + fe Zo) A. 


Solving eqns. (17)—(22) for fo, 8, 1 81 and fo, Z. one gets 





fo 2 fA (fo — 9° — 1 24) 
Ao _ 2 fo (fo — 2) Bulld (fo — 0" — UI (25) 
4h. = — thy (fo — 2) + bi No (fo — 9 — 1) (26) 
FS = thy + (fo — 2) hall fo — 0 — 1 ee) 
ass = —[hs (fo — A) + hulllgo (fo — 4? — .. (28) 
G82 — tha + ha (fo — Nil fo — 0° — 1 voe(29) 


Equations (24)—(29) are six simultaneous differential equations info, 80, Si, 81, fz and 
g,. Boundary conditions for fo, 80; fi, 81, fy and g, are obtained from relations (5)— 
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(6), i. e. the jump conditions across the shock front. At the shock where A = 1, we 
have 


...(30 
8 =1 g, = 0 8, = 0. (30) 


We have integrated the differential equations (24) — (29) subject to the boundary 
conditions (30), using Runge-Kutta method of fourth order the results are shown in 
the Figs. 1 to 4. 
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Fic. 1. Variation of f versus } for weak shocks. 


4. DISCUSSION AND THE CONCLUSIONS 


In Figs. | and 2, we have shown the variation of the parameters f and g versus 
A for AP = 0.001, 0.005 and 0.01 respectively. This is the case of weak shocks. In 
Table I, values of shock pressure P2 are shown for various values of AP. 


It is seen from Fig. 1 that particle velocity ratios u 


/U =f (A) increases continuously 
as A decreases from | to 0.1, 


value of f being approximately zero at A = 1 for weak 
shocks. In Fig. 2, variation of density ratio g (A) is shown. Itis seen that for Ap 
= 0.001, g first increases and when A — 0.19, it starts decreasing. The trend of the 
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Fic. 2. Variation of g versus 4 for weak shocks. 
TABLE I 

ee ——————EeEeEeEeEeEeEeEeEEEEEEEE—— 

Weak Shock Strong shock 

Ap p2(Kb) Bo p2(Kb) 

» = Se oe a 
0.001 0.0243 0.150 5.6062 
0.005 0.1229 0.250 12.6925 
0.010 0.2492 0.350 24.6650 


EE 


density variations in the present problem is similar to that of piston problem‘. In Fig. 
3 fis plotted versus A for the case of strong shock waves. ForA = 1, value of f = u./ 
U = (8 — 1)/5, which first decreases, then starts increasing exponentially, as A decre- 
ases from | to zero. Similar trend is found in the variation of g for strong shocks. 


Once conditions at the shock front are known, variation of fluid parameters is 
known from shock front to the point of explosion. Conditions at the shock front are 


functions of shock radius, which is given by Singh et ai.*. 


In the present work we have ignored the presence of gas bubble at the centre. 
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Fic. 3. Variation of f versus » for strong shocks. 


ACKNOWLEDGEMENT 


Authors are thankful to Professor P. L. Sachdev, Indian Institute of Science, for 
helpful discussions and to Dr SM Sharma, Scientist E, for the encouragement. Thanks 
are also due to Dr R Natarajan, Director, CASSA, for giving permission to publish 
this paper. We are also grateful to referee for his useful comments; which helped us to 
bring this paper in the present form. 


REFERENCES 


1. N.N. Kochina and N. S. Melnikova, J. Appl. Math. Mech. 23 (1959) 123. 
N. N. Kochina and N. S. Melnikova, Proc. Stokloy. Inst. Math. 87 (1966), p. 31. 


. 


a ae 


a 


SHOCK WAVES IN WATER 201 


A. K. Madan et al., Establishment of Acquarium Technique; TBRL Report No. 239/83 (Re- 
Stricted), 1983. 


M. P. Ranga Rao and B. V. Ramana, Int. J. Engng Sci. 11 (1973) 1317. 
P. Singh, Bull. Cal. Math. Soc. 63 (1971), 87-96. 

V. P. Singh and M.S. Bola, Indian J. pure Appl. Math. 7 (1976), 1405. 
V. P. Singh, Indian J. pure appl. Math. 7 (1976), 147. 


V. P. Singh, A. K. Madan, H.R. Saneja and Dal Chand, Proc. Indian Acad. Sci. 3 (1980), 
169. 


125 


075 
ao 


050 


025 





025 0.50 0.75 10 
A —o 


Fic. 4. Variation of g versus» for strong shocks. 


APPENDIX A 
Derivation of the Equation (10) 


Substituting equation (9) in (8) we get 


(1+ Ap) Ap? _—  K_ ..(A.1) 
(1 + A(1 — 5) R 
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where 


K = 34QJ R® [47 p, a’. .. (A.2) 


Taking logs and differentiating (A.1) with respect to R we get 


al Sif Oe — b) + Ap? (1 — b)] [1 + 2 AP 
ok = ap Lt AP (2 > b) ened TN ites 
i as Ap]. (A.3) 


Expanding the right handside and rearranging the coefficient of AP and Ap’ etc. 
we get 


Ap BAP 


oR R 
Where 


B = By + B, AP + B, Ap? + ... 
Bo = 3/2 
B, = —3 (2b — 1)/4 
Py = 3 (4b — 1)/8. 
2. Derivation of Equation (13) and (14) 


I 


Since entropy variations are negligible in underwater shocks, we have 


ees ( 2) oP nn OP 
or op }, 


——— SE pe 


or : or 
where 

c* = @ p/éPp), 
is the sound velocity in compressed water. 
Using equation (5) we get 


2. 2 _ a(S +50 —1) 


ep, [3 —6 (8 — Dp 





In equations (1) and (2) independent parameters r and ¢ are transformed to non-di- 
mensional parameter A and Ap by the following operators 


Cas oe Ve ee 
ot Raa R oAp 
0 7) 


l 
Or: ee eee 
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Using (A.5) alongwith (A.4), (9) (10) and (12) in equations (1) and (2) we get equa- 
tion (13)—(14) after some simplifications. 


3. Comparision with Ref. Ranga and Ramana‘ 
Although problem in reference 4 is dealt using similarity methods, but trend in the 
variations of fluid parameters behind the shock in piston problem and explosion pro- 


blem should be comparables. The similarity parameter A in ref.‘ is same asA of our 
paper, as follows. In reference‘ equation (33) and (18) are as 


A = (A/P,)8? rt-8 ...(33) 


ro = at, ...(18) 
Eliminating ¢ from (18) and (33) we get 

A = (A/P1)8!? & (r/rs) 
Now at the shock front r = r,, A = | 


ra (A/p1)5!? aan |, 
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